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Abstract 

DELPHI data collected at centre-of-mass energies up to 208 GeV have been 
analysed to search for charginos, neutralinos and sfermions in the framework of 
the Minimal Supersymmetric Standard Model (MSSM) with R-parity conserva- 
tion. No evidence for a signal was found in any of the channels. The results of 
each search were used to derive limits on production cross-sections and particle 
masses. In addition, the combined result of all searches excludes regions in the 
parameter space of the constrained MSSM, leading to limits on the mass of the 
Lightest Supersymmetric Particle and other supersymmetric particles. 
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1 Introduction 

Supersymmetry (SUSY) [1] is at present one of the most attractive possible exten- 
sions of the Standard Model (SM) and its signatures could be observed at LEP through 
a large variety of different channels. This paper presents searches for the pair-production 
of charginos, neutralinos, sleptons and squarks. The searches were performed and inter- 
preted in the most model-independent way possible in terms of production cross-sections 
and masses. The results were interpreted in the framework of constrained SUSY models, 
with the different search channels complementing each other in constraining the param- 
eter space. 

The data collected by the DELPHI experiment in e+e~ coUisions at centre-of-mass 
energies {^/s) up to 208 GeV were used. No signal was found in any of the channels, and 
limits were set at 95% confidence level (CL). 

The paper is organised as follows. In section 2 the basic supersymmetry framework is 
described: the general phenomenology is discussed in section 2.1, and implications of more 
constrained models used for interpreting the data are given in section 2.2. The DELPHI 
detector is described in section 3, and in section 4 the data sets and event generators are 
reviewed. In section 5 the general analysis framework is described, the analysis methods 
are briefly mentioned and the specific searches for sleptons, squarks, charginos and neu- 
tralinos are discussed. The results of each search are separately presented and interpreted 
in section 6. In section 7 the results are combined and interpreted in the framework of 
constrained SUSY scenarios with gravity-induced breaking of supersymmetry (SUGRA). 
A brief summary is given in section 8. 

Previous results published by DELPHI can be found in references [2] to [9]. 

2 SUSY framework 

The searches presented in this paper were performed in the framework of the Minimal 
Supersymmetric extension of the Standard Model (MSSM) [1]. i?-parity ^ conservation 
is assumed, implying that the Lightest Supersymmetric Particle (LSP) is stable and 
SUSY particles ( "sparticles" ) are pair-produced. In addition, they decay directly or 
indirectly into the LSP. In this paper the lightest neutrahno, Xi, is assumed to be the 
weakly interacting LSP, which escapes detection giving signatures of missing energy and 
momentum. 

The searches for sparticle production were developed with minimal assumptions, and 
the selections employed depended primarily on the masses of the particles involved. In 
particular, the sensitivity of the searches depends on the visible energy released in the 
decay process. In direct decays into the LSP, this visible energy is largely determined 
by the mass difference (AM) between the decaying sparticle and the LSP. In indirect 
(cascade) decays, other mass differences can also be important. 

The MSSM has a large number of free parameters in addition to the SM ones. The 
most model-independent interpretation of the results is in terms of the masses and cross- 
sections explicitly involved, for each production channel. A common interpretation of 
the results from the various searches can also be performed and used to exclude regions 
of the model parameter space. This, however, requires a manageable number of free 
parameters. For this reason, assumptions must be made and specific scenarios defined 

^R-parity is a multiplicative quantum number defined as R = (— where B, L and S are the baryon number, 
the lepton number and the spin of the particle, respectively. SM particles have R=+l while their SUSY partners have 
R= -1. 
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for such an interpretation. The general phenomenofogy of the searches will be discussed 
in section 2.1, followed by a description of the implications of more specific scenarios in 
section 2.2. 

2.1 General phenomenology 

Squarks and Sleptons 

The "sfermions", squarks and sleptons, are the scalar partners of the SM fermions. 
The left- and right-handed chiral states of each SM fermion, fL and fa, have as SUSY 
partners two scalars, usually labelled fL and fR. 

Sleptons and squarks could be pair- produced at LEP via e+e~ annihilation into Z/7, 

leading to f^fR or fLfL final states. Selectrons could also be produced through t-channel 
neutralino exchange. The selectron cross-section depends critically on the neutralino 
mass, and destructive interference can make it very small. The t-channel contribution 
also introduces the possibility of clCr and crCl production. 

If the unification of sfermion masses at a high mass scale typical of Grand Unified 
Theories (GUT) is assumed, smaller masses and cross-sections for a given universal mass 
parameter arc typically expected for the partners of right-handed fermions. Under this 
assumption, the kinematic accessibility of first and second family sfermions at LEP de- 
pends only on their assumed common mass at the unification scale. Squarks are in general 
expected to be heavier than sleptons. However, for sfermions of the third family the large 
Yukawa couplings lower the masses, as large mixing between left and right states may oc- 
cur. In this case, the lighter mass states of third family sleptons and squarks, fi (stau), bi 
(sbottom) and ti(stop), are candidates for the lightest charged supersymmetric particle. 

In large regions of the SUSY parameter space the dominant decay of the sfermions 
is to the corresponding fermion and the lightest neutralino, f ^ f^^- In the case of the 
stop, the decay t —>■ txi is not kinematically allowed at LEP, and the dominant two-body 
decay channel is expected to be t ^ cxI (t hxf being disfavoured by existing limits 
on the chargino mass). If mo < M^^, the three-body decay ti hii> may compete with 
the cXi decay. 

Thus final state topologies with a pair of acoplanar ^ leptons or jets and missing 
energy are the relevant ones in the search for sleptons and squarks, respectively, and the 
total energy of the detectable final state particles (and thus the sensitivity of the search) 
depends primarily on the mass difference between the sfermion and the LSP. 

Gharginos and Neutrafinos 

In the MSSM there are four neutralinos, = 1,4 (numbered in order of increasing 
mass) and two charginos x^, i = 1, 2 which are linear combinations of the SUSY partners 
of neutral and charged gauge and Higgs bosons (gauginos and higgsinos). The fightest 
states can be mainly gaugino or higgsino, or strongly mixed (for similar gaugino and 
higgsino mass parameters). 

Neutralinos (charginos) could be pair-produced at LEP via s-channel Z (Z/7) ex- 
change or i-channel exchange of a selectron (sneutrino). The i-channel contribution can 
be important if the slepton is light. The interference with the s-channel diagram is 
constructive in the case of neutralinos, but destructive in the case of charginos. 

'^In this context, acoplanar means that the direction of one of tfie Icptons/jets is not in tfie plane defined by the direction 
of the other lepton/jct and the beam hnc. Throughout this paper, the acoplanarity angle between two particles or jets is 
defined as the complementary (with respect to 180°) of the angle between their directions in the plane transverse to the 
beam direction. 
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Expected decays are x° — > x?f f and xf ~^ Xi^ f'- If th^ sfermions are heavy, these 
decays proceed via Z or W exchange. However, sfermion emission may also contribute 
if the sfermions are hght, increasing the partial width for decays into the corresponding 
fermions, and two-body decays into f f can dominate if kinematically allowed. The one- 
loop decay X2 ~^ Xil can be important in specific regions of the parameter space when 
other decays are suppressed. 

Thus the sfermion mass spectrum may significantly affect both the production cross- 
section and the decay modes of charginos and neutralinos, and many final state topologies 
are possible. 

In the case of chargino pair-production, the final state is expected to be four jets if 
both charginos decay hadronically, two jets and one lepton if one chargino decays into 
ii^Xi: ^-nd leptons only if both charginos decay into leptons. The branching ratio of 
xf ^ X2^ f' can be sizable, in particular in the regions of the parameter space where 
X2 X1I is important. In this case, the above topologies are accompanied by photons. 

If the mass difference AM between the chargino and the LSP is very small the visible 
energy released in the decay is very small, making the signal hard to detect. The simulta- 
neous production of a photon by initial state radiation (ISR) can be used to explore such 
regions, as this allows a very efficient background rejection (at the expense of a small 
signal cross-section). Still smaller mass differences imply a long lifetime of the chargino, 
which can then be identified as a heavy stable charged particle or one with a displaced 
decay vertex. 

In the case of the detectable XiXt neutralino production channels {i.e. excluding XiXi)^ 
the most important signatures are expected to be acoplanar pairs of jets or leptons with 
high missing energy and momentum. Although X1X2 X1X3 expected to dominate 
in most of the parameter space, channels like X2X3 X2X4 must also be considered for 
a complete coverage. These give rise to cascade decays with multiple jets or leptons in 
the final state, possibly accompanied by photons from X2 ~^ Xi7- 

2.2 Constraining the parameter space 

In this paper the results are combined and interpreted in the framework of constrained 
SUSY scenarios with gravity-induced SUSY breaking. To make the model more predic- 
tive, the unification of some parameters at a high mass scale typical of Grand Unified 
Theories (GUT) can be assumed [10]. The MSSM parameters and the assumptions that 
can be relevant in the interpretation of the results are listed below: 

• tan /3, the ratio of the vacuum expectation values of the two Higgs doublets; 

• /Li, the Higgsino mixing mass parameter; 

• Ml, M2, Ms, the f/(l) X SU{2) x SU{3) gaugino masses at the electroweak (EW) 
scale: when gaugino mass unification at the GUT scale is assumed, with a common 
gaugino mass mi/2, the relation between Mi and M2 is Mi = |M2 tan^ 9w ~ O.5M2, 
with the assumption that Mi/ai {i — 1,2,3) is renormalization group invariant, 
fulfilled at the one-loop level; 

• Mf, the sfermion masses: under the assumption of sfermion mass unification, mo is 
the common sfermion mass at the GUT scale; 

• the pseudoscalar Higgs mass, tua, on which the masses in the Higgs sector depend: 
if scalar mass unification is assumed, rriA at the EW scale can be derived from mo; 

• the trilinear couplings, Af, determining the mixing in the sfermion families: the 
third family trilinear couplings are the most relevant ones, A^-, A^,, At, and under 
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the assumption of universal parameters at high mass scale there is a common trilinear 
coupling A, to which A^, and At at the EW scale can be related. 

Mass mixing terms at the EW scale given by mr^AT — /xtan/3), mjy^Aij — /lianP) and 
mt{At — /x/tan/3) are considered for f, b and t, respectively. The mass splitting grows 
with the mixing terms, and for large |//| this can give light fi and bi states if tan/3 is 
large, or a light ti for small tan /9. 

In the model referred to as the "Constrained MSSM" (CMSSM) in the following, 
sfermion and gaugino mass unification are assumed. The parameter set is then reduced 
to M2, mo, tan (3, n, rriA and three Af couplings. This is the model considered in section 7. 

Further tightening the assumptions, the individual Af couplings can be replaced by a 
universal coupling A and rriA can be related to the other parameters by assuming scalar 
(inchiding Higgs bosons) mass unification. Requiring in addition the correct reproduction 
of the EW symmetry scale, which fixes the absolute value of n, defines the minimal MSSM 
with gravity- induced SUSY breaking (mSUGRA). 

The direct results of the searches are first derived and presented in the most model- 
independent way possible. Under the assumptions described above, the results are then 
used to constrain the SUSY parameters. Presently, the strongest constraints on SUSY 
models come from the MSSM Higgs searches [1 1] . 

Chargino production is the most important direct SUSY detection channel for large 
regions of the parameter space. However, if sfermions are light (corresponding to a small 
mo scenario), or if the parameters take particular values, the chargino production cross- 
section can be greatly suppressed or undetectable final states can dominate (in particular 
for small mass sphttings) . The most relevant of these cases are the following: 

Large mo 

For large mo, the sfermions are heavy and have little infiuence on the observable 
phenomenology. The chargino pair-production cross-section is large and the chargino is 
excluded nearly up to the kinematic limit. The search for charginos in non- degenerate 
scenarios has been applied down to AM = —M^o = 3 GeV/c^. The region AM < 
3 GeV/c^ is covered by the search requiring an ISR photon and by the searches for long- 
lived charginos. Moreover, small values of AM occurring in scenarios without gaugino 
mass unification [12] have also been investigated. 

At small tan/3 (tan/3 < 1.2), neutralino searches can extend the sensitivity of the 
searches beyond the kinematic limit for chargino production. This concerns the region 
of (small) negative n and M2>60 GeV/c^. In particular, searches for neutralino cascade 
decays are crucial for investigating the M2 < 120 GeV/c^ region, where the X1X2 cross- 
section is small. 

Small mo 

If mo is small, light sfermions affect the chargino and neutralino production cross- 
sections. In particular, the decrease of mo causes the chargino production cross-section 
to drop in the region where the gaugino components dominate (small M2 and large |/x|). 
Down to mo — 200 GeV/c^, the cross-section remains large enough to allow chargino 
exclusion nearly up to the kinematic limit. For smaller mo, the neutralino production 
cross-section is very much enhanced, and neutralino searches become sensitive instead. 

If mo is very small (mo ~ 100 GeV/c^), sleptons can be sufficiently fight to affect 
drastically the decay patterns of charginos and neutralinos, and nearly invisible final 
states can become dominant in some cases. However, for such small mo and small M2 
(below 200 GeV/c^) sleptons can also be searched for in direct pair-production. 
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If M-± > rriy and the mass difference M-± — mj> is small, the chargino decay chain 

xf — > i>£ — s> z^Xi^ is dominant, and leads to an experimentally undetectable final state 
(the only visible final state lepton has very low momentum). However, in this case the 
search for selectrons can be used to constrain the sneutrino mass (under the assumption 
of unification) and thus the chargino mass. 

It can also happen, in scenarios with large mixing among sfermions and large tan/3 and 
M2, that "blind spots" occur in the chargino detection sensitivity due to xf — >fi/— >rxii^ 
with a small mass difference Mf—M^o. In this case, X1X2 or X2X2 production with xl^^'^ 
are the only detectable channels. A specific search was designed for this case. 

3 Detector description 

The DELPHI detector is described in detail in [13]. The central tracking system 
consisted of a Time Projection Chamber (TPC), supplemented by a system of silicon 
tracking detectors and drift chambers. These included the Vertex Detector (VD), closest 
to the beam pipe, the Inner Detector (ID) and the Outer Detector (OD). These were 
situated inside a solenoidal magnetic field of 1.2 T, parallel to the beam axis. The 
average momentum resolution for charged particles in hadronic final states was in the 
range Ap/p^ ~ 0.001 - 0.01 (GeV/c)-^. 

The electromagnetic (EM) calorimeters were symmetric around the plane perpendicu- 
lar to the beam {6 = 90°) ^, with the High density Projection Chamber (HPC) extending 
from 88.7° to 43.1° (barrel region), the Forward Electromagnetic Calorimeter (FEMC) 
from 36° down to 9°, overlapping with the Small angle Tile Calorimeter (STIC), the 
DELPHI luminometer, which covered the range 1.7° < 9 < 10.6°. The region of poor 
electromagnetic calorimetry at a polar angle close to 40° was instrumented by scintillators 
(hermeticity taggers) [14] to tag photons. 

The Hadron Calorimeter (HCAL) covered 98% of the solid angle. Muons with mo- 
menta above 2 GeV/c penetrated the HCAL and were recorded in sets of Muon drift 
chambers located in the barrel (MUB), forward (MUF) and surround ^ (MUS) regions of 
the detector. 

4 Data samples and event generators 

During the year 2000 DELPHI collected data in the centre-of-mass energy range from 
201.5 to 208.8 GeV. The average centre-of-mass energy was < -\/i>~ 206 GeV and the 
total integrated luminosity amoTinted to about 224 pb""^. In 1999 (1998) about 227 pb~^ 
(158 pb~^) were collected at centre-of-mass energies around 192, 196, 200 and 202 GeV 
(189 GeV). 

The data collected in the years 1999-2000 were analysed in the searches presented in 
this paper. In some cases the 1998 data were re-analysed for consistency with the im- 
proved methods now presented, as detailed in the description of the analyses. Combina- 
tion with earlier results was performed to obtain cross-section upper limits and excluded 
regions in the model parameter space. Details for each channel are given in section 6. 

^In DELPHI, a right-handed Cartesian coordinate system is used with the z direction defined by tiic direction of the 
electron beam, and the x-axis pointing towards the centre of the LEP ring. The origin is at the centre of the detector. The 
polar and azimuthal angles 6 and (j) are defined with respect to the z axis and <^ = corresponds to the a;-direction, and 
the coordinate r is defined in the usual way as r = ya;^ + j/^ . In this paper, polar angle ranges are always assumed to be 
symmetric with respect to the 6 = 90° axis. 

^The region between the barrel and end-cap parts of DELPHI not covered by the MUF and the MUB chambers. 
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Quality requirements on the status of the subdetectors most relevant for each analysis 
were applied, generically based on the status of the main tracking devices. Muon cham- 
bers or calorimeters were considered in channels where muon or electromagnetic shower 
detection was crucial. The luminosity loss was at most of the order of a few percent and 
was taken into account in the analyses. 

On September 2000, sector 6 of the TPC (corresponding to 1/12 of the TPC 
acceptance) failed beyond repair. This required modifications of the pattern recognition, 
and affected the quality of charged track reconstruction. Thus special care had to be 
taken for each search when analysing the data collected without TPC sector 6. The 
accumulated integrated luminosities with and without a working sector 6 are respectively 
164 pb~^ and 60 pb~^at average centre-of-mass energies around 206 GeV. 

In order to increase the sensitivity for a discovery, the data collected in 2000 were 
divided into 4 regions of centre-of-mass energy as given in table 1. 



^/s region of analysis 




< > (GeV) 


1 


Sector 6 on 


v/i< 205.75 GcV 


78.3 


204.9 


2 


Sector 6 on 


205.75 <yi< 207.5 GeV 


78.8 


206.7 


3 


Sector 6 on 


207.5 GeV< ^ 


7.2 


208.0 


4 


Sector 6 off 


all ^/s 


60.0 


206.5(*) 



Table 1: Definition of the ^/s regions used to analyse the data collected in 2000. The 
three energy bins for the period in which the detector was fully operational are referred 
to in the text by their approximate average ccntrc-of-mass energies: 205, 207 and 208 
GeV. The "Sector 6 off" data are referred to as 206.5 (*). 

To evaluate the signal efficiencies and background contaminations, simulated events 
were generated using several different programs. 

The background process e+e~ — > qq(n7) was generated with PYTHIA 6. 125 [15]. For 
A«"^A«~(7) and t+t"(7), DYMU3 [16] and KORALZ 4.2 [17] were used, respectively. The 
BHWIDE generator [18] was used for Bhabha events. Simulation of four-fermion final 
states was performed using EXCALIBUR [19] and grc4f [20]. 

Two-photon interactions giving hadronic final states were generated using TWOGAM [21], 
PHQJET [22] and PYTHIA 6.143 [15], while leptonic final states were generated using the 
generator of [23], including radiative corrections for the e"^e~ /z^//"" and e'^'e" r+r" final 
states. 

SUSYGEN 2 . 2004 [24] was used to generate chargino, neutralino, slcpton, and sbottom 
signal events and to calculate cross-sections and branching ratios for these channels. For 
the nearly mass-degenerate case, the chargino decays were modelled with the results of 
the computations of reference [12]. Stop events were generated according to the expected 
differential cross-sections using the BASES and SPRING program packages and taking spe- 
cial care in the modelling of the stop hadronisation [25]. 

In all cases except for stop generation, JETSET 7.4 [15], tuned to LEPl data [26], was 
used for quark fragmentation. 

The generated signal and background events were passed through the detailed simula- 
tion of the detector, DELSIM [13], and then processed with the same reconstruction and 
analysis programs as the real data. The faster DELPHI simulation code SGV [27] was 
also used for signal efficiency studies in some analyses. SGV is a model of the detector 
response which has been tuned to the data independently from DELSIM. For charged 
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particles the perigee parameters resulting from the physical process generation were mod- 
ified to take into account parametrized effects of multiple scattering, detector resolutions 
and acceptance. For neutral particles geometric acceptance and parametrized calorimet- 
ric resolutions were taken into account. Bremsstrahlung and photon conversion in the 
tracking system were also simulated. 

The numbers of simulated events from different background processes were several 
times the numbers in the real data. 

5 Descriptions of the analyses 

The analyses described below can be divided into two stages. The first stage was very 
similar for all searches and consisted of the selection of charged and neutral particles 
followed by an event preselection. These are described in sections 5.1.1 and 5.1.2. In 
the subsequent stage, the analyses differed according to the characteristics of the vari- 
ous signals. This is described in sections 5.2 to 5.5. In most searches several different 
topologies were considered, accounting for the different possible final states. Particle 
identification and reconstruction algorithms common to several searches are described in 
section 5.1.3. Different analysis techniques were chosen for the various searches: some 
analyses were based on successive requirements on individual event variables ("sequen- 
tial cut analyses"), others used multidimensional techniques based on likelihood ratios or 
neural networks. These techniques are briefly discussed in section 5.1.4. 

As discussed above, the sensitivity of the searches depends on the mass difference 
between the produced sparticles and the LSP, which determines the visible energy released 
in the process. Typically, for small AM the signatures of sparticle production arc similar 
to those of two-photon interactions (77 events). For high AM sparticle production events 
resemble four-fermion final states such as W'''W~ and ZZ. For intermediate AM values, 
the background is composed of several SM processes (in particular two-fermion ones). 
The coverage of all the relevant AM regions often requires the combination of different 
searches or the optimisation of the selection criteria separately in each AM interval. 
Sparticle searches for very low values of AM are particularly challenging and required 
different preselections. 

5.1 Basic selections and techniques 
5.1.1 Particle selection 

The following quality requirements were applied to the charged and neutral particles 
observed in the detector. 

Charged particles were required to have momentum p above 100 MeV/c and below 
0.75^/s, a relative momentum error less than 100%, and to extrapolate back to within 5 
cm of the main vertex in the transverse (r0) plane and 10 cm/sin in the longitudinal 
[z) direction. Similar but more stringent criteria were applied to particles whose tracks 
extrapolated to the TPC, but which gave no signal in the TPC Whenever tracks recon- 
structed with TPC information were required in the subsequent analyses, at least five of 
the 16 pad rows had to contribute hits. 

In the stau analysis (see section 5.2.3) charged particle tracks were required to have 
TPC information, or all three of the detectors VD, ID and OD used in the reconstruction 
of the track. In addition, only tracks with polar angle 9 > 15° were kept. 
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In the nearly mass-degenerate chargino search (see section 5.4.2) there was no lower 
bound on the momentum for tracks at polar angle above 25°, while p > 150 MeV/c was 
required otherwise. In addition, different impact parameter requirements were applied (6 
cm and 12 cm in the transverse and logitudinal directions, respectively). 

Energy clusters in the calorimeters were taken as neutral particles if not associated to 
a charged particle and if above an energy threshold which was 900 MeV for deposits in the 
hadron calorimeter and ranged from 300 to 500 MeV for deposits in the electromagnetic 
calorimeters (depending on the region of the detector). Cuts removing clusters created 
by radioactivity in the lead of the HPC or by particles from cosmic ray showers were also 
applied. 

5.1.2 General event preselection 

In the general preselection, events were kept if there were at least two charged particles, 
at least one of them had a transverse momentum above 1.5 GeV/c, and the transverse 
energy of the event ^ exceeded 4 GeV. This rejected mostly two-photon interactions (for 
which most of the energy is deposited in the forward regions of the detector), zero- or 
one-prong final states (like e+e" —>■ jj, e^y —>■ ej) and beam-related backgrounds (such 
as beam-gas interactions). 

For chargino searches in nearly mass-degenerate scenarios a different preselection was 
used: at least two charged particles were required, as well as one isolated electromagnetic 
cluster with transverse energy above a y^-dependent threshold close to 5 GcV and a 
mass recoiling against it above 90 GeV/c^. Two-photon and beam-gas backgrounds were 
reduced by rejecting events with a large fraction of the detected energy in the forward 
region of the detector. 

5.1.3 Pcirticle identification and reconstruction algorithms 

The following criteria for particle and event classification were common to the different 
searches. 

Particle jets were reconstructed using three different approaches: 

• The DURHAM [28] algorithm was used to cluster the particles into a fixed number of 
jets: two or four. 

• The LUCLUS [15] algorithm was applied with the critical distance set to d^om — 
10 GeV/c or djoin — 2.5 GeV/c. The final number of jets is, thus, variable (and 
lower in the first case). 

• A specific algorithm optimised for the low multiplicity jets resulting from r decays 
was used for f pair-production searches. This method considered all possible ways 
of clustering the charged particles in the event into groups, always requiring the 
invariant mass to be below 2 GcV/c^. Clearly identified leptons were considered as 
a single group, except for pairs of oppositely charged, well identified electrons close 
together which were allowed to be grouped with other particles, since they could 
come from a converted photon. 

If possible, the event was clustered into two groups with invariant mass below 

2 GcV/c^. If no such combination existed, the one with the smallest number of 
groups was kept. When more than one way of obtaining two groups both with in- 
variant mass below 2 GeV/c^ was found, the grouping yielding the lowest sum of 
masses was retained. Once the best grouping of the charged tracks was found, it 

^Thc transverse energy is defined as the sum of the absolute values of the transverse momenta of all particles in the 

event: c^iPTi- 
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was attempted to associate the neutrals in the event to the particle groups. Also 
in this step, the grouping yielding the lowest sum of masses was chosen. However, 
as very few long-lived neutral hadrons are expected in r-decays, neutral hadronic 
clusters were not included in the groups, but treated as isolated neutrals. Also elec- 
tromagnetic clusters which could not be merged into any of the groups without the 
invariant mass exceeding 2 GeV/c^, were left as isolated neutrals. In addition, a 
special procedure was applied to identify and correct for neutral clusters that were 
likely to be either bremsstrahlung photons or a shower induced by an electron that 
was not correctly assigned to the track by the reconstruction program. 
The charge of each cluster was taken as the sum of the charges of all its particles. 

The thrust variable used in several analyses throughout this paper was computed using 
all the particles meeting the particle selection requirements. 

Isolated leptons or photons in the event are often very important in distinguishing 
signal and background. In the present searches, the isolation criteria depended on the 
multiplicity of the event: 

• In low multiplicity searches, charged particles were classified as isolated if the total 
charged energy, excluding the energy of the particle itself, within 10° of the track 
direction was below 2 GcV. Slightly tighter cuts in the impact parameters (1 cm 
and 5 cm in rcj) and z, respectively) and in the momentum error were also applied. 

• In high multiplicity searches, a photon was considered isolated if its angular separa- 
tion from any neutral or charged particle was greater than 15°. A lepton was tagged 
as isolated if its angular separation from all the jets (computed without the lepton 
using the LUCLUS algorithm with djom = 40 GeV/c) was greater than 20°. 

The identification of a track as a muon, electron, or hadron was "tight" , "loose" , or 
"veto" (or none at all). Tight identification was unambiguous. A particle could simulta- 
neously be loosely identified as several different species. Excluded particle species were 
vetoed. 

The identification of muons was provided primarily by the DELPHI standard algorithm 
described in [13], which relies on the association of charged particles to signals in the 
muon chambers and the HCAL. 

Electron and photon identification was performed by the algorithm described in [29] 
which combines deposits in the EM calorimeters with tracking information and takes 
possible radiation and interaction effects into account by a clustering procedure in an 
angular region around the main shower. In the f analysis (see section 5.2.3), the clustering 
procedure was not used for tracks in the barrel region of the detector, since it tended to 
treat charged pions from r decays wrongly, if they were accompanied closely by neutral 
pions. In low multiplicity topologies a very loose electron identification based on the 
ratio E /p between the energy deposited in the EM calorimeter and the momentum of the 
associated charged track was also used. 

Below, whenever the identification level is not specified, it is imphed to be "tight" for 
electrons and "loose" for muons. In the case of electrons, "tight" identification basically 
adds some isolation requirements to the identification ones. 

e/7 separation inside the acceptance of the STIC luminometer was performed on a 
statistical basis, using the veto information of the two planes of the scintillator counters 
placed in front of it. 

In the stau analysis a particle was considered as a tight hadron if it was not classified 
as a muon or electron and had an associated energy in the hadron calorimeter exceeding 
50% of its momentum, or else was considered as a loose hadron if it had hadronic energy 
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associated and it was not tightly identified as a muon or electron. If both the electro- 
magnetic and hadronic energies were small (less than 1 and 6 GeV, respectively), and 
the difference between the hadronic energy and the track momentum was above 10 GeV, 
the particle was assigned both the loose muon and loose hadron code. 

Decays of b-quarks were tagged using a probabilistic method based on the impact 
parameter of tracks with respect to the main vertex [30]. stands for the corresponding 
probability estimator for tracks with positive impact parameter, the sign of the impact 
parameter being defined by the jet direction. The combined probability Vcomb included 
additional contributions from properties of reconstructed secondary vertices. 

All searches made use of the information from the hermeticity taggers [14] to reject 
events with photons in the otherwise insensitive region at polar angles around 40°. If 
there were active taggers not associated to reconstructed jets, the event was rejected if the 
tagger was located in the direction of the missing momentum. In the stau analysis, where 
neutrinos from tau decays made the estimation of the direction of the missing momentum 
unreliable, events containing active and isolated taggers were rejected irrespective of the 
direction of the missing momentum. 

5.1.4 Analysis techniques 

Likelihood ratio method 

In the likelihood ratio method used, several discriminating variables arc combined 
into one on the basis of their one-dimensional probability density functions (pdf's). If 
the variables used are independent, this gives the best possible background suppression 
for a given signal efficiency [31]. For a set of variables {xj}, the pdf's of these variables 
are estimated by normalised frequency distributions for the signal and the background 
samples. We denote the pdf's of these variables ff{xi) for the signal events and f,P{xi) 
for the background events submitted to the same selection criteria. The likelihood ratio 

" fS(x) 

function is defined as C-jz — j| fB{x) ' ^^^^^^ with jC-j^ > jC-jz^^^ are selected as candidate 

signal events. The choice of variables and the value of C-ji^^^^ were optimised using 
simulated event samples by minimising the signal cross-section that was expected to be 
excluded at 95% confidence level in the absence of a signal. 

Neural networks 

A neural network provides a different way of defining one discriminating variable from 
multidimensional distributions of event variables given as inputs. The neural network 
used below (see section 5.3) contains three layers of nodes: the input layer where each 
neuron corresponds to a discriminating variable, the hidden layer, and the output layer 
which is the response of the neural network. The program used in the squark analysis 
was SNNS [32]. A "feed- forward" architecture is implemented and the "back-propagation" 
algorithm is used to train the network with simulated events. An independent validation 
sample was also used not to overtrain the network. A way of enhancing the efficiency 
of the network without increasing too much the number of its parameters is to define 
a separate output node for each type of event that the neural network should separate. 
More details are given in section 5.3. 
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5.2 Slept on searches 

Supersymmetric partners of electrons, muons and tans were searched for. In this 
paper, data collected at centre-of-mass energies between 189 and 208 GeV were analysed, 
and were combined with previous results [2]. 

The track selection and the general event preselection described in sections 5.1.1 and 
5.1.2, respectively, were used. The analyses were then performed in two stages. Firstly, 
a loose selection was used to obtain a sample of low multiplicity events. Events with less 
than ten charged particles and a visible invariant mass above 4.5 GeV/c^ were retained 
for further analysis. Different selections were then applied in each of the three channels. 

5.2.1 Selectron seEirches 

To search for selectrons, the general topology required was two acoplanar electrons 
and missing energy. All candidates with exactly two well reconstructed and isolated 
particles (according to the definition described in section 5.1.3), oppositely charged and 
with momentum above 1 GeV/c, were first selected. One of the two charged particles was 
required to be tightly identified as an electron, and the event was rejected if the other 
was identified as a muon. 

At this stage of the selection the sample consisted mainly of Bhabha and two-photon 
events. Satisfactory agreement was observed between the data and simulated background, 
as shown in figure 1 ^. 

A series of tighter cuts reduced the SM background further. As two-photon events 
are predominantly at low polar angles and with low momentum, the visible energy was 
required to be greater than 15 GeV, the energy deposited in the low angle STIC calorime- 
ter less than 4 GeV and the total transverse momentum with respect to the beam axis 
greater than 5 GeV/c. To reduce the number of Bhabha events an upper limit on the 
visible energy of 100 GeV was imposed, while also requiring that the neutral energy not 
associated to the charged particle tracks be less than 30 GeV. Events were also rejected 
if there were more than four neutral clusters in total, each with energy above 500 MeV. 
Bhabha events are coplanar and with a large opening angle, hence it was demanded that 
the opening angle between the two tracks be lower than 165° and the acoplanarity be 
greater than 15°. Constraints were also imposed on the momenta of the two particles, 
requiring that both tracks had momentum greater than 2 GeV/c. 

5.2.2 Smuon searches 

Smuon pair-production with decays to muon plus neutralino is expected to give acopla- 
nar muons and missing energy. All candidates with exactly two well reconstructed and 
isolated particles (according to the definition given in section 5.1.3), oppositely charged 
and with momentum above 1 GeV/c, were first selected. At least one of the particles 
had to be loosely identified as a muon. It was further required that neither particle be 
identified as an electron. The selected sample consisted mainly of two-photon events and 
fair agreement between real data and simulated background was observed (see figure 2; 
the slight discrepancies visible in figure 2(d), in the tail of the two-photon distribution 
at about 40 GeV/c, and for high momenta, are in kinematic regions rejected in the next 
step of the analysis). 

To reduce the SM background further, a series of tighter cuts were applied. To remove 
two-photon events, the visible energy was required to be greater than 10 GeV. Also, the 

®In order to show the diflferent background contributions, the largely dominant Bhabha background was suppressed in 
these plots by demanding that the opening angle between the two tra<;ks was below 176° . 
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energy in the STIC had to be less than 1 GeV. It was also demanded that the transverse 
momentum be greater than 5 GeV/c. To remove e+e" -^ii^ ^~ events, an upper limit of 
120 GeV on the visible energy was imposed, whilst also requiring the unassociated neutral 
energy to be less than 10 GcV, with no more than two neutral clusters. This background 
was further suppressed by accepting only events in which the opening angle between the 
tracks was less than 165° and the acoplanarity was greater than 15°. To reduce t-channel 
W pair contamination, events were rejected if the positively charged muon was within 
40° of the e"*" beam direction, or the negatively charged muon was within 40° degrees of 
the e~ beam direction. 

5.2.3 Stau searches 

Events with two acoplanar taus and high missing energy can be the signature of 
stau pair-production. Due to the scalar nature of the stau, the two taus are produced 
centrally in the detector. To select this topology, the particles in the events were grouped 
into clusters according to the algorithm described in section 5.1.3. Events with exactly 
two particle clusters (possibly accompanied by isolated neutral particles) were considered 
further if they contained at least one charged particle with momentum above 1 GeV/c 
and a relative error less than 30%. This particle had to be isolated (no neutral or 
charged particles in a cone of 20° half-angle around it), or above 30° in polar angle, 
and its calorimetric energy could not exceed the momentum by more than three times 
the expected error on the calorimetric measurement. A comparison of data and simulated 
SM background at this stage is shown in figure 3. It was further required that there 
were no more than six charged particles in the event and that the total charge was or 
±1. 

Tight and loose electrons, muons and hadrons were defined as described in sec- 
tion 5.1.3. In addition, if a particle had a loose identification for one species, and was 
vetoed for the other two, it was considered as tightly identified. If a track pointed to a 
gap in the electromagnetic calorimetry, it was considered as a loose electron. Tracks with 
no identification information were treated as loose electrons. Neutral clusters passing the 
criteria of section 5.1.1 were used if their angle to the beam was above 15°. Furthermore, 
no identified hadronic secondary interactions inside the tracking system were allowed. 

Beyond this point, the analysis differed depending on whether a stau with mass above 
or below mz/2 was searched for. 

Search for staus with large mass 



To suppress the two-photon background, it was required that the total transverse mo- 
mentum imbalance exceeded 4 GeV/c, the total calorimetric energy below 30° in polar 
angle did not exceed 20% of the beam momentum and the total momentum of the event 
was within the region 9{T,p) > 30°. 

To reduce the background from radiative return to the Z, none of the clusters was 
allowed to have a total momentum {p'^^'^) above 70% of the beam momentum, the mo- 
mentum of isolated photons had to be less than 10% of the beam momentum, and the 
acoplanarity was required to be above 12°. To reject e+e~ — >• Z/7 — >• r+r^ events where 
the decay of one r yielded visible products with large momentum, while the decay of the 

^In order to show the contribution of all classes of background, the largely dominant Bliabha background was suppressed 
in these plots by demanding that the missing transverse momentum was above 4 GeV/c, and that the acoplanarity was 
above 1°. 
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other r yielded soft products, the momentum transverse to the thrust axis was required 
to exceed 0.7 GeV/c. 

This selection was supplemented by cuts that depended on the region of the (Mf ,Mlsp) 
plane considered. At any given point in the plane, it was demanded that the ac- 
cepted events in that point were kinematically compatible with the corresponding sig- 
nal: the maximal momentum the r can have in the lab frame, neglecting m,-, is 

Pmax = ^ (l - [Mi^sp/Mff) (l + a/1 - AMf'^/s^, which is also the end-point of the 

spectrum of the visible r - jet momentum. Hence p^^^ (the larger of the two jet mo- 
menta) was required to be less than Pmax- 

For large AM= — Mlsp, the remaining background from two-photon events can be 
removed by requiring large transverse momentum imbalance. At smaller values of AM 
such a cut would greatly reduce the signal detection efficiency. Therefore, if AM was 
below (above) 20 GeV/c^ the total missing transverse momentum (/^ ) had to exceed 0.8 
(1.2) times the maximum transverse momentum a 77 event could have without one of 
the beam-remnant electrons being deflected into the STIC {i.e. by an angle greater than 
0max = 1.82°). This limit depended on the centre-of-mass energy: j;!^"^ = ^/s sm9„iax/ii + 
sinOraax) = 0.031 y/s ^. In addition, in the region with AM below 20 GeV/c^ there should 
be no calorimetric energy below 30°, and O^Ep) had to exceed 45°. 

At this stage of the analysis, the dominant remaining background at large AM were 
W^W~ events. Only events with both W bosons decaying leptonically were still present. 
Electrons and muons in W'''W~ events might come either directly from the W-decay, 
or indirectly from r-decays. In the former case, which is dominant, the momentum of 
the detected lepton tends to be higher than in the signal, where all electrons and muons 
would be indirect. In the latter case, the momentum spectrum of background and signal 
events is similar. 

In order to reduce the background from W'''W~ events with leptons from direct W 
decays, it was demanded that the highest momentum of any tightly or loosely identified 
lepton in the event was less than (0.1 AM -|- 0.6)P^^j^, where Pwmin is the lowest 
momentum a lepton in the decay W iu can obtain in the lab-frame if the W is on- 
shell. In the region of large AM, it was also demanded that there be no more than one 
tightly identified electron or muon in the event. 

To suppress further the W^W" background, and in particular the component with 
leptons from indirect decays, the events were analysed as if they were indeed W'''W~ 
events. The 9 angle of the W yielding the cluster with positive charge was estimated 
{9w+) and selections in the (^vK+5Pfe^^)-pla'ne were used to discriminate the signal from 
this background ^. The 9w+ was estimated as follows. If neither W decayed to a tau, 
9w+ can be calculated exactly (albeit usually with two- fold ambiguity). The W decay to 
a muon or an electron could be distinguished from the decay to the r by the presence of 
a single track cluster with at least a loose lepton identification and a momentum above 
the lowest possible momentum for the charged lepton from the decay of an on-shell W. 
Decays that did not fulfil these requirements were assumed to be decays to r. If one of 
the W decayed to a r the solution for 9\y+ is approximate, since the momentum of the 
T is unknown. In this case, the momentum of the initial r was estimated as the average 
momentum of taus from W decay, calculated for tau momenta above the measured jet 

**Thc missing transverse momentum was estimated in four different ways: from the transverse momentum of the two jets, 
as that of all particles passing the quality cuts, as that of all reconstructed particles except those identified as bremsstrahlung 
photons, and as that of all reconstructed particles. The cut was applied on the smallest of these. This gave stability against 
possible errors in reconstruction and against the presence of noise or cosmics. 

^At this stage of the analysis about 99 % of the events had a cluster with charge equal to +1 (this selection was not 
aplied to events where there was no cluster with positive charge). 
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momentum. The direction of the momentum of the tau was taken as the measured jet 
direction. 

The W pair- production process tends to have higher p^^^ and 9w+- In addition, p'j^^fi 
and 9w+ are correlated in such processes, while they are independent in f production. 
Thus the signal selection criteria in the -plane were: 

C -3.0(x- 0.325)+^ if x< 0.325 
0W+ < \ A if 0.325 <x < 0.52 

[ -2.1(x-0.52)+^ if x> 0.52 

where 9w+ is in radians and x — pj^^h/ Pbeam, and ^4 is a constant chosen to be 1.6 (2.1) for 

AM below (above) 20 GeV/c^. The boundary of the selected region in the {dw+:P'M^h)~ 
plane thus defined closely follows a curve of constant ratio between the probability density 
functions for signal and background (the likelihood ratio). 

Search for light staus without coupling to the Z 



A light stau can be exchidcd to a large extent using LEPl results, as discussed in sec- 
tion 6.3. This is however not possible when the coupling to the Z vanishes (the stau 
mixing angle gives the minimum cross-section) . The large mass analysis described in the 
previous section looses its efficiency for stau masses below 15 GeV/c^. This is mainly due 
to the fact that the staus are highly boosted at such small masses, failing the acoplanarity 
cut. 

Therefore a specific search was required for small Mf at the minimal cross-section 
mixing angle. Two search regions were identified: one optimised for very small masses, 
Mf below 10 GeV/c^, referred to as the "very small mass analysis", and one for larger 
masses, optimized for Mf above 10 GeV/c^and AM between rrir and 4 GeV/c^ referred 
to as the "small mass analysis" (the large mass analysis shows good sensitivity for AM 
above 4 GeV/c^ down to Mf = 15 GeV/c^). 

The signal events in the relevant kinematic regions are characterised by containing two 
taus at large angles to the beam, and being softer than two-fermion events but slightly 
harder than 77 events. Due to the sizable boost of the staus, the two jets tend to be 
rather back-to-back. Two-tau events were selected as described in the previous section, 
with the additional requirement that the topology was either 1-prong and 3-prong, or two 
1-prongs. In the latter case, there should not be two tightly identified leptons of the same 
fiavour in the event. To select central events, the polar angle of the most energetic particle 
in each detector hemisphere had to be above 50°, and the sine of the polar angle of each 
jet had to be above 0.8. The acoUinearity was required to be above 0.4°, p^^/pbeam had 
to be below 90%, and the total reconstructed mass of the visible system had to be above 
4.5 GeV/c^. There had to be no energy in a 30° cone around the beam axis. 

In the small mass region, 0{llp) had to be above 55° in polar angle, the visible mass 
had to be below [15(AM - m^) - {Mf - 25)/3 - (AM - m^)(Mf - 25)/2 + 15] GeV/c^ 
and ft had to exceed max(0.05, (AM- m^)/4) p^™. 

In the very small mass region, it was required that ft was above O.OlpJf™, pinlr/Pheam 
was above 15%, and 9{'Lp) was above 15° to the beam. The acoUinearity had to be below 
15°. 
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5.3 Squark searches 

Supersymmetric partners of top and bottom quarks were searched for. The data 
collected at centre-of-mass energies from 189 to 208 GeV were analysed. The dominant 
decays of the stop and sbottom squarks are assumed to be ti — > cxi and bi — > hxi, 
respectively, and the final topology is two acoplanar jets and missing energy. In the 
non-degenerate scenario (AM > 10 GeV/c^), the neural network analysis has already 
been presented in [3]. This analysis has been extended down to AM = 5 GeV/c^. In 
addition, a new analysis based on a sequential cut approach has been developed to search 
for stops nearly degenerate in mass with the LSP, investigating AM values between 2 and 
10 GcV/c^. Moreover, this analysis has been extended to AM values up to 20 GeV/c^ in 
order to cross-check the non-dcgcncratc analysis. 

The track selection and the general event preselection described in sections 5.1.1 and 
5.1.2, respectively, were used. 

5.3.1 Non-degenerate scencirios 

To select hadronic events, the number of charged particles reconstructed with TPC 
information was required to be greater than three, and the energy in the STIC to be less 
than 70% of the detected energy. The polar angle of the thrust axis had to be above 20°. 
The following event quality cuts were then applied. The percentage of good tracks, the 
ratio of the number of charged particle tracks after the particle selection to the number 
before, had to be greater than 35%. In addition, the scalar sum of charged particle 
momenta reconstructed with TPC information was required to be greater than 55% of 
the total energy in the event, and the total number of charged particles to be greater 
than six. 

To remove radiative return events, the energy of the most energetic neutral particle 
was required to be less than 40 GeV. Additional cuts were then applied to restrict the 
selection to events with missing energy. The transverse missing momentum had to be 
greater than 4 GeV/c, the polar angle of the missing momentum had to be above 20° 
and the energy in a 40° cone around the z axis was required to be less than 40% of the 
total detected event energy. Finally, the visible mass of the events was required to be less 
than 95 GeV/c^. 

The number of events selected by this preselection was 2178 for 2143±8 expected 
(combined data from ^/s =189 to 208 GeV). Figure 4 shows a comparison between 
data and simulated events. At this level, for AM > 10 GeV/c^, stop signal efficiencies 
ranged from 20% to 70% depending on the mass difference between the stop and the 
neutralino. Sbottom efficiencies were quite similar except at low AM where, for example, 
the efficiency for Mj^^=90 GeV/c^ and M^o=85 GeV/c^ was close to zero, because the b 
quarks are produced almost at rest. 

The final selection of events was performed using neural network techniques (see section 
5.1.4). Separate searches were made for two different ranges of AM: AM > 20 GeV/c^ 
and 5 < AM < 20 GeV/c^. Events were forced into two jets using the Durham algorithm. 
The neural network structure was as follows. There were ten input nodes (variables), ten 
hidden nodes (in one layer) and three output nodes. The ten input variables were: the 
ratio between the transverse missing momentum and the visible energy, the transverse 
energy, the visible mass, the softness defined as MjetiZ-Ejeti + Mjet2/-Ejct2) the acollinear- 

ity, the quadratic sum of the transverse momenta of the jets (Pf^^)'^ + (P^^*^)^, the 
acoplanarity, the sum of the first and third Fox- Wolfram moments, the polar angle of 
the missing momentum and finally the combined b-tagging event probability. For each 
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AM window a neural network with three ouput nodes was trained to discriminate the 
signal from the combined two-fermion and four-fermion backgrounds, and from the 77 
interactions leading to hadronic final states. 

Although the three output nodes proved useful in training the network, the selection 
was made according to the output of the signal node only. Figure 5 shows the number of 
events as a function of the signal efficiency for the two mass analysis windows of the stop 
and the sbottom searches. The number of events in the data is in agreement with the SM 
background predictions over the full range of neural network outputs. The optimisation 
of the final cuts was performed by minimising the confidence level of the signal hypothesis 
expected in the absence of a signal. [33] . 

5.3.2 Nearly mass degenerate scenarios 

Due to the large Yukawa coupling (see section 2.1), the stop (ti) and the sbottom (61) 
can be fight and nearly degenerate in mass with the LSP. The effective coupfing of the 
stop to charm and neutralino results from loops and is thus small. In addition, the width 
of the decay t cXi is proportional to Mi(l — M^t)^ / Mi^^ , and therefore proportional 
to AM. So if AM gets small enough, the stop acquires a sizable lifetime and may form 
a quasi-stable (decaying inside the tracking volume) or even stable stop hadron (see [34] 
for this case). The current analysis focusses on a stop decaying promptly into a charm 
particle and the LSP. 

The event preselection required, in addition to the criteria described in section 5.1.2, 
that not more than 30% of the total visible energy was carried by particles with tracks 
seen in the VD and ID only. To eliminate Bhabhas and leptonic 77 backgrounds, the 
charged multiplicity was required to be greater than five. The 77 background was further 
suppressed by requiring the energy in a forward cone of 30° around the beam direction to 
be at most 40% of the total visible energy and smaller than 2 GeV, and that no energy 
was deposited in the STIC calorimeter. To avoid the relatively low hadronic energy 
region, where the 77 background is not well reproduced by the MC, the total transverse 
charged energy was required to be greater than 7 GeV, the total transverse energy of 
tracks reconstructed with TPC information and the total transverse momentum had to 
be greater than 4 GeV and 3.5 GeV/c, respectively, and the number of tracks with TPC 
information had to be at least four. 

The agreement between data and MC after this preselection is shown in fig- 
ure 6(a) to 6(d). Figure 6 also demonstrates that the two-fermion and four-fermion 
backgrounds dominate at this stage of the selection. The smaller contribution from two- 
photon interactions can be reduced at this level using for example the different polar 
angle distribution with respect to the signal. 

A further selection was performed in order to reduce the remaining backgrounds. 
Events having mainly barrel activity were selected. This was achieved by requiring that 
the energy within a cone of 60° around the beam direction was less than 10 GeV and that 
the polar angle of the missing momentum was above 45°. Most of the remaining two- and 
four-fermion background was rejected by demanding that the transverse momentum of 
the most energetic particle was less than 10 GeV/c and that the total transverse energy 
was less than 40 GeV. Finally, the total transverse momentum was required to be greater 
than 5 GeV/c and the scaled acoplanarity ^° greater than 20°. This cut removed most of 
the remaining background from two-photon processes. The agreement between data and 
simulation after this selection is shown in figure 6(e,f). 

'^^The scaled acoplanarity is the acoplanarity of the two jets multiplied by the sine of the minimum angle between a jet 
and the beam axis. 
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5.4 Chargino searches 

5.4.1 Non-degenerate scenEirios 

The search for charginos in the non- degenerate scenarios covers the case when the mass 
difference AM= Af-±— M^o is above 3 GeV/c^. In order to take all possible signatures of 
chargino decays into account, events were divided into four mutually exclusive topologies: 

• the ii topology, with no more than five charged particles and no isolated photons; 

• the jji topology, with more than five charged particles, at least one isolated lepton 
and no isolated photons; 

• the jets topology, with more than five charged particles and no isolated photons or 

leptons; 

• the rad topology, with at least one isolated photon. 

The signal events selected in a given topology are mostly events from the corresponding 
decay channel, but events from other channels may also contribute. For instance, for low 
AM (and thus low visible energy) some events with hadronic decays are selected in the 
ii topology, and some mixed decay events with the isolated lepton unidentified enter into 
the jets topology. This migration effect tends to disappear as AM increases. This effect 
was taken into account in the final efficiency and limit computations. 

The signal events were simulated using 132 combinations of xt ^-^id Xi masses for 
nine chargino mass values (M-± ^ 103, 102, 100, 98, 94, 85, 70, 50 and 45 GeV/c^) and 
with AM ranging from 3 GeV/c^ to 80 GeV/c^. A total of 264000 chargino events (2000 
per mass combination) was generated. The kinematic observables (acoplanarity^^, i?vis! 
ft J etc.) of the signal events were studied in terms of their mean value and standard 
deviation, and six AM regions were defined, each containing signal events with similar 
properties (see table 2). 

In each of these 24 windows (four topologies, six AM regions), a likelihood ratio 
function {C-r, see section 5.1.4) was defined. The variables {xi} used to build the C-ji 
functions in the present analysis were [4]: the visible energy (£^vis), visible mass (Myis), 
missing transverse momentum {ft), polar angle of the missing momentum, number of 
charged particles, total number of particles, acoplanarity, acoUinearity, ratio of electro- 
magnetic energy to total detected energy, percentage of total energy within 30° of the 
beam axis, kinematic information concerning the isolated photons and leptons and the 
two most energetic charged particles, and finally the jet characteristics. 

The generation of these 24 likehhood ratio functions was performed as follows: 

• The signal distributions of all the variables {xi} were built with signal events gen- 
erated with parameter sets giving rise to charginos and neutralinos with masses in 
the corresponding AM region. For each AM region the events were classified ac- 
cording to the above topological cuts. The background distributions were built with 
background events passing the same topological cuts. 

• Preselection cuts [35], different for each AM region, were apphed in order to reduce 
the backgrounds with largest cross-section (two-photon interactions and Bhabha 
events) and to generate the pdf's. The total background was reduced to 5% of 
the one passing the general event preselection (section 5.1.2). The pdf's were then 
generated as mentioned in section 5.1.4. Figure 7(8) shows the distributions of some 
event variables for the jets, ii and rad topologies for the 2000 data and simulation 
with the TPC sector 6 on (off). 

^^To compute the acoplanarity and acoUinearity the particles were forced into two jets by the DURHAM algorithm. 
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• To reduce statistical fluctuations, a smoothing was performed by passing the 24 sets 
of pdf 's for signal and background through a triangular filter [36] . 

• In each window all the combinations of the pdf 's were tested, starting from a minimal 
set of four variables. Every combination defined an C-ji function and an H-ricuri ^ 
described in section 5.1.4, using the single channel formula [37]. The parameters 
entering this computation were the number of expected background events and the 
window efficiency of the chargino selection, defined as the mean efficiency of the 
chargino-neutralino mass sets belonging to the investigated window Figures 7(d) 
and 8(d) show the good agreement obtained between real and simulated events as a 
function of the likelihood ratio cut, for 25 < AM < 35 GeV/c^ in the jj£ topology. 

• The combination of variables corresponding to the lowest excluded cross-section 
defined the L-ji function and the C-Rcur window. 



AM regions 


1 


3<AM< 5 GeV/c' 


2 


5<AM< 10 GeV/c' 


3 


10<AM< 25 GeV/c^ 


4 


25<AM< 35 GeV/c^ 


5 


35<AM< 50 GcV/c" 


6 


50 GcV/c^<AM 



Table 2: Definition of the AM (mass difference between the chargino and lightest neu- 
tralino) regions for the chargino search in non-degenerate scenarios. 



5.4.2 Nearly mass-degenerate scenarios 

The search for charginos in the nearly mass-degenerate scenarios uses several different 
techniques, depending on the lifetime of the chargino. which in turn depends on the mass 
difference AM between the chargino and the lightest neutralino (this is the only relevant 
dependence, at least in the heavy slepton hypothesis). When AM is below the mass of 
the pion, the chargino lifetime is usually long enough to let it pass through the entire 
detector before decaying. This range of AM can be covered by the search for long-lived 
heavy charged particles. For AM of a few hundred MeV/c^ the chargino can decay inside 
the main tracking devices. Therefore, a search for secondary vertices or kinks can be used 
to cover this region. As the mass difference increases, the mean lifetime shortens until 
the position of the X\ decay can hardly be distinguished from the main event vertex. In 
this case, the tagging of an energetic ISR photon can help in exploring the AM region 
between a few hundred MeV/c^ and 3 GeV/c^. The selection criteria are similar to the 
ones used in the analysis of previous data, which have been described in [6] . 

Search for long lived charginos 

Long lived charginos can either be "quasi-stable" (decay outside the tracking system) 
or decay "visibly" inside the tracking devices. 

The search for heavy stable charged particles is described in [5] . The method used to 
identify heavy stable particles relied on the ionisation loss measurements in the TPC and 

^^The efficiency of one chargino-neutralino mass set is defined as the number of events satisfying Cfi > C-jz^ut divided 
by the total number of chargino events satisfying the topological cuts. 
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on the absence of signal in the DELPHI Cherenkov radiation detectors (RICH). Heavy 
stable charged particles crossing the detector would be seen in the tracking system and 
have as distinctive signature the absence of Cherenkov radiation and an anomalous energy 
loss in the TPC. Three different search windows were used in the search for heavy stable 
charginos: 

• the charged particle had momentum above 15 GeV/c, and no photons in any of the 
two radiators of the RICH (liquid, refractive index n — 1.28, and gas, n — 1.0015) 
were associated to the track; 

• the charged particle had momentum above 5 GeV/c, high ionisation loss in the TPC, 
and no signal in the gas RICH; 

• the charged particle had momentum above 15 GeV/c, a TPC ionisation loss not 
exceeding 70% of the expectation for a proton, and no signal in the gas RICH. 

A fourth search window considered in [5] was not included, in order to treat the two 
hemispheres of the event independently. 

If a heavy charged particle decays inside the central tracking devices (at a radius 
between 10 cm and 1 m) then both the incoming and the outgoing track can be recon- 
structed, and the angle between the tracks can be calculated. Such a search for kinks was 
originally designed to search for long-lived staus in the Gauge Mediated SUSY Breaking 
scenario [5]. A similar technique was applied to search for mass-degenerate charginos, 
with some specific features needed because the visible decay products carry very little 
momentum in the nearly mass-degenerate case. Details of the selection criteria can be 
found in [6] . Here only a brief and qualitative summary of the most important selection 
cuts is given. 

A set of rather loose general requirements was applied in order to suppress the low 
energy background (beam-gas, beam- wall, etc), two-photon, e"'"e~ and hadronic events. 
For each event passing the preselection cuts, all the charged particles were grouped in 
clusters according to their measured point closest to the interaction vertex. A cluster with 
only one track with momentum above 20 GeV/c was considered as a possible chargino 
candidate if it was compatible with a particle coming from the interaction point. For 
each single track cluster fulfilling the above conditions, a search was made for a second 
cluster possibly formed by the decay products of the and defining a secondary vertex 
or kink with the chargino candidate. 

Reconstructed secondary vertices could also be the result of particles interacting in 
the detector material, or having a particle trajectory reconstructed in two separate track 
segments. Additional requirements rejected these backgrounds in the events with an 
acceptable secondary vertex [5]. Finally, for an event to be accepted, at least one charged 
particle had to be found in each hemisphere (defined by the plane containing the beam 
spot and perpendicular to the line connecting the beam spot to the kink). 

The search for events with tracks at large impact parameter described in [5] was 
not possible in this case: events with only two extremely soft charged particles with 
large impact parameter are difficult both to trigger on and to discriminate from machine 
related noise. Such events were however considered if a high pt ISR photon was present, 
as explained in the next section. 

Search for charginos with ISR photons 

The visible particles resulting from the decay of a chargino nearly mass-degenerate 
with the LSP have typically little energy and momentum. The trigger efficiency is low 
for such events, and there is a very large background from two-photon events. The ISR 
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photon tag improves detectability and, if the transverse energy of the photon is above 
a threshold which depends on the minimal polar angle acceptance of the experiment, it 
rejects most of the two-photon background. 

After the preselection, which was summarised in 5.1.2, the following requirements were 
applied to the data and simulation samples. 

• There had to be at least two and at most six good charged particles passing the 
quality criteria (see section 5.1.1), and no more than ten tracks in total. 

• The transverse energy of the ISR photon candidate was required to be greater than 
(E^)""'^ ~ 0.03 • v^. 

• The mass recoiling against the photon had to be above 2M-± — 5M, where the term 

Al 

SM takes into account the energy resolution in the electromagnetic calorimeters. 

• The photon had to be isolated by at least 30° with respect to any other charged or 
neutral particle in the event. 

• The sum of the energies of the particles with polar angles within 30° of the beam 
axis (£"30) was required to be less than 25% of the total visible energy. If the photon 
itself was below 30°, it was the ratio (£^30 — E^)/{Eyis — Ej) that was required to be 
below 0.25. 

• If the ISR photon candidate was detected in the very forward calorimeter STIC, it 
must not be correlated with a signal in the scintillators placed in front of the STIC. 

• (^vis ~ Ey) / ^/s had to be below a kinematic threshold which depended on AM and 
on M-± (and in any case below 6%). 

• The ratio of the absolute value of the missing transverse momentum over the total 
transverse energy had to be above 0.40/c if AM > 300 MeV/c^, and above 0.75/c 
for smaller AM. 

• If AM > 1 GcV/c^ , at least two charged particles in the event had to be consistent 
with coming from the beam interaction region. 

Distributions of some of the variables used in the final selection are shown in figure 9 for 
data, simulated SM background, and simulated signal events. Although there is a certain 
overall qualitative agreement of the various distributions, there is already an excess of 
data. On the other hand, the two-photon generators used in the simulation lack the 
events which have small 77 invariant mass, and in some cases (namely, 77 — > e"^e~) no 
ISR generation is implemented at all. Moreover, background processes such as beam-gas 
interactions arc not included in the simulation. As in previous publications [6] , the most 
likely explanation of such disagreement is therefore a deficit of simulated background 
events rather than an excess of data from possible new physics. As no attempt will be 
made in the following to account for the backgrounds missing in the simulation, the limits 
that will be obtained are conservative. 

5.5 Neutralino searches 

The neutralino searches were designed to cover both xtXi production with xt 
Xi+ff, with a signature of acoplanar jets or leptons, and channels of the type xtx^j with 
k or i > 2, which can lead to neutralino cascade decays. To maximise the sensitivity 
several searches were used, covering different topologies, namely: 

• a search for acoplanar jet events, as from X1X2 with X2 ~^ Xi'i^i 

• a search for acoplanar lepton events, as from with X2 ~^ Xi^'^^" or 

X2 x?/^^/^~; 
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• a search for multijet events, as from XiXji'^ = l)2,j = 3,4 with Xj X2CL<1 and 

decaying to x?qq or x57; 

• a search for multilepton events for the corresponding decays to lepton pairs; 

• a search for cascade decays with tau leptons, e.g. X2X1 production with X2 ~^ 

and f x?^; 

• a search for double cascade decays with tau leptons, e.g. X2X2 production with the 
same X2 decay chain as above. 

The different searches, briefly described below, were designed to be mutually exclusive 
in order to allow easy combination of the results. Thus events selected in the likelihood- 
based searches for acoplanar leptons or jets of section 5.5.1 were explicitly rejected in the 
searches described in the subsequent sections. The track selection and the general event 
preselection described in sections 5.1.1 and 5.1.2, respectively, were used. The LUCLUS 
algorithm with djom = 10 GeV/c was used for jet clustering in the analyses described 
below, with the exceptions of the likelihood ratio acoplanar jets/leptons and double tau 
cascade searches, as explicitly mentioned in the corresponding sections (see section 5.1.3 
for details on the jet clustering algorithms). 

Data collected at centre-of-mass energies between 192 and 208 GeV were analysed, 
and were combined with previous results [7]. In the search for staus in X2X2 production, 
all data recorded at centre-of-mass energies between 189 and 208 GeV were analysed. 

5.5.1 Acoplancir jets and acoplancir leptons seeirches 

As mentioned above, the acoplanar jets and acoplanar leptons topologies are dominant 

in most of the parameter space. For these search based on the likelihood ratio 

method was performed and the sequential cut analyses described in [7] were used as a 
cross-check. 

The characteristics of the neutralino decays are mainly determined by the value of 

AM, here defined as the mass difference between the heavier of the produced neutralinos 
and the LSP (x?)- The total energy of the visible final state particles, E^is, was used 
to distinguish between regions of different signal and background characteristics in the 
optimisation of the selections. 

Likelihood Ratio analysis 

The first step of the analysis was to preselect the events dividing them into three 
mutually exclusive topologies: ee, mi, and qq. The ee and nn topologies were defined as 
having exactly two isolated lepton candidates (see section 5.1.3). In the ee topology, at 
least one of these had to be a tightly identified electron, and neither identified as a muon. 
Similarly, the events in the /i/j, topology were required to contain at least one isolated 
loose muon candidate and no isolated electron. The qq topology was defined as events 
with more than five charged particles and no isolated photons or leptons. 

In the second step, aimed at removing the dominant SM background processes, events 
which fulfilled all of the following criteria were selected: 

• the polar angles of the most energetic neutral and charged particles were required 

to be above 10°; 

• the missing transverse momentum had to exceed 2 GeV/c, or 4 GeV/c if the visible 
energy was less than 30 GeV; 

• both the acoplanarity and acollinearity had to be greater than 3°; 

• the total visible energy had to be lower than 0.75-^. 



22 



The first two requirements remove the bulk of the 77 events and off-momentum beam 
electrons. The third and fourth reject two-fermion processes. The last requirement 
removes mainly four-fermion events. 

In the qq topology an additional selection was applied to suppress further the large 
qq(7) background. This was based on the jets reconstructed using the Durham algorithm 
and forcing the number of jets to two. If the invariant mass of the event was within 
40 GeV/c^ of the Z mass, the acoplanarity of the two jets was required to be at least 10°. 

For different values of visible energy E^ig (typical of different values of AM), the 
kinematic properties of the signal were studied in terms of the mean value and standard 
deviation of several event variables. Five £^vis regions were defined, each containing signal 
events with similar properties and SM background composition. These regions are given 
in table 3, together with the corresponding dominant SM background. 



-Evis regions 


Main SM l)kg. 


1 


5<^vis< 20 GcV/c^ 


77 


2 


20<Evis< 50 GcV/c' 


77, 2-fermions 


3 


50<E^is< 70 GeV/c' 


2-fermions 


4 


70<E,i,< 110 GeV/c" 


2- and 4-fermions 


5 


110 GeV/c^<£;vis 


4-fermions 



Table 3: Definition of the visible energy regions of the neutralino search in the acoplanar 
leptons and jets topologies and corresponding dominant SM backgrounds. 

In the last step of the analysis, for each of the 15 windows (three topologies, five ii^vis 
regions) thus defined, a likehhood ratio function was computed. The variables used in 
the hkelihood definition are listed below: 

• global variables (all topologies): visible energy, transverse energy, missing momen- 
tum, energy and direction of the most energetic charged and neutral particles, trans- 
verse momentum with respect to the thrust axis, polar angle of the missing momen- 
tum, thrust value, thrust direction and acoplanarity; 

• variables specific to the qq topology: jet directions, energies, widths, invariant mass 
and ^cut value; 

• variable specific to the ee and /i/i topologies, namely the invariant mass of the two 
charged particles. 

For each variable used in the likelihood ratio function, the one-dimensional probability 
density function was defined according to the procedure described in sections 5.1.4 and 
5.4.1, and events with C-ji > C-ji^^^ were selected as candidate signal events. 

Some event variable distributions for real and simulated data before the likelihood 
selection are shown in figure 10. A fair agreement between data and the SM expectation 
is found. A low visible energy event is present in the ee topology (see figure 10(d)), in a 
region where low but non-zero background was expected from the SM (about 0.5 events 
with visible energy below 30 GeV expected from four-fermion processes). No significant 
disagreement between data and MC is found in any of the topologies. 

Seciuential analysis 

The detailed selection criteria for the selection based on sequential cuts are given in 
reference [7] and have not been changed. This analysis was used as a cross-check of the 
likelihood ratio results. At the final selection level, criteria optimised for different AM 
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regions were designed. They were used as independent selections in the derivation of the 
results. 

5.5.2 Multijet seeirch 

The multijet search was optimised for cascade decays of neutralinos with large mass 
splittings, giving high energy jets. Events with energetic photons, characteristic of the 
decay X2 ~^ Xi7) were treated separately. Events with a photon signature were selected 
on the basis of reconstructed photons with a polar angle above 20°, isolated by more than 
20° from the nearest charged particle. If there was only one such photon its energy was 
required to be between 10 GeV and 40 GeV; if more than one photon was present, at 
least two had to have energy greater than 10 GeV. 

This selection was similar to the acoplanar jet selection of the sequential analysis, but 
required a rather large transverse energy and allowed any number of reconstructed jets. 
Events selected by the searches for acoplanar leptons or acoplanar jets were explicitly 
rejected. 

The detailed selection criteria are similar to the ones described in [7]. Background 
studies based on data of 1998 and an improved energy-flow reconstruction motivated 
several changes in the selection procedure. Three selection stages (preselection, interme- 
diate selection and final selection) are defined and used in the figures and tables. In the 
following, the most important steps of this search are summarised. 

At preselection level, at least five charged particles (passing the track selection de- 
scribed in section 5.1.1) were required, at least one of them with a transverse momentum 
exceeding 2.5 GeV/c. The transverse energy of the event had to be greater than 25 GeV, 
the visible energy was required to be less than 0.65^/s, and the missing momentum had 
to be less than OAy^/c. There were several requirements aimed at selecting events 
with jets which were not dominated by single particles with large reconstructed energy. 
Figures ll(a,b) and 12(a,b) show the distributions of the visible mass divided by the 
centre-of-mass energy for real data and simulated background events passing the above 
selection. 

At the intermediate selection level, radiative return to the Z, two-photon, and Bhabha 
background was reduced by excluding events with a neutral particle whose energy ex- 
ceeded 60 GeV or with more than 40% of the visible energy within 30° of the beam 
direction. Also, the transverse momentum had to exceed 6 GeV/c. The total momentum 
and the most energetic shower in the event were both required not to be close to the 
beam direction. A comparison of the cpt/\/s distributions for real data and simulated 
background following the above selection is shown in figures ll(c,d) and 12(c,d). The ex- 
cess in data visible in 11(c) and 12(c) was found to be consistent with (low acoplanarity) 
radiative return to the Z events with misreconstructed missing momentum. In 12(c) the 
excess is larger. The reason for this is that in the year 2000 the missing momentum 
reconstruction suffered from the problems with sector 6 of the TPC and from the fact 
that the distortions in the TPC data had not yet been fully corrected for at the time of 
this analysis. The phase space region in question was removed by criteria imposed later 
in the selection, and the excess has no impact on the analysis results. 

At the final level of the selection the acoUinearity and scaled acoplanarity (in a forced 
2-jet configuration) had to be greater than 30° and 10°, respectively. To reject W^W~ 
background it was required that there be no charged particle with a momentum above 
30 GeV/c, and no isolated lepton above 10 GeV/c or above 4 GeV/c with an isolation 
angle greater than 20°. 



24 



For events with a photon signature, the mass recoihng against the visible system was 
required to exceed 20% of the centre-of-mass energy and the scalar sum of momenta 
reconstructed with TPC information had to be less than 60% of the visible energy. 

For the complementary sample, without a photon signature, several of the criteria 
above were made stricter in order to reject Z7 events. The recoil mass was required to 
exceed 40% of \/i and the energy deposited in electromagnetic calorimeters had to be 
less than 40 GeV. No isolated neutrals with energy greater than 20 GeV were allowed 
and the average momentum of particles with tracks reconstructed in the TPC had to be 
less than 4 GeV/c. 

Figures ll(e,f) and 12(e,f) show the distribution of the scaled acoplanarity for real 
and simulated data after final selections, while figures 12(g,h) show cpt/\/s for data with 
TPC sector 6 on. A satisfactory agreement between data and SM simulation is found. 
Four events are observed while three were expected from SM simulation. 

5.5.3 Multilepton search 

The multilepton search is sensitive to cascade decays involving leptons, which can 
dominate if there are light sleptons. This search was described in [7] and is briefly 
summarised here. 

At the preselection level, well reconstructed low multiplicity events with missing en- 
ergy and missing mass were selected. In particular, the total visible energy including 
badly reconstructed tracks (not passing the track selection described in section 5.1.1) 
was required to be less than 140 GeV, the number of charged particles was required to 
be at least two and at most eight, and events with more than four neutral particles were 
rejected. Figures 13(a) and 14(a) show a comparison between the visible mass divided 
by the ccntrc-of-mass energy for real and simulated events passing the preselection. 

The selection at the intermediate level served mainly to reject Z7, two-photon, and 
Bhabha events by requiring significant transverse momentum and transverse energy (pt> 
8 GeV/c, -E'^s> 25 GeV). The distributions of cp-^lsfs for real and simulated data, fol- 
lowing the intermediate selection are compared in figures 13(c) and 14(c). 

At the final selection level events with two or more charged particles were subjected to 
criteria designed to reject Bhabha events. In addition, the charge asymmetry for the two 
most energetic such particles was used to reject W pairs decaying leptonically. Events 
with four or more tracks were clustered into jets, and those with exactly two jets were 
rejected if their scaled acoplanarity was less than 15°. Figures 13(e) and 14(e) show the 
acoplanarity distribution for real and simulated data at this level. 

5.5.4 Asymmetric tau cascade search 

The tau cascade search is sensitive to x5x2 production with X2 ~^ ^'^ ^ ~^ Xi''") 
where the second r produced has very low energy. This search was described in [7] and 
is briefiy summarised here. 

At the preselection level, well reconstructed low multiplicity events with missing energy 
and missing mass were selected. The selection was the same as for the multilepton search 
(section 5.5.3), with the additional requirement of no more than two reconstructed jets. 
At least two of the charged particles had also to satisfy stricter criteria on reconstruction 
and impact parameters. The distributions of the visible mass divided by the centre-of- 
mass energy for real and simulated data at this level are shown in figures 13(b) and 
14(b). 
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At the intermediate selection level, the highest and second highest momenta of charged 
particles were required to be below 40 GeV/c and 25 GeV/c, respectively, and at least one 
charged particle had to have a transverse momentum above 2.5 GeV/c. The criteria to 
reject Z7, two-photon, and Bhabha events were similar to those used in the multilepton 
search, except for the removal of the requirement on the transverse energy. Figures 13(d) 
and 14(d) show the distributions of cp'^j \fs for real and simulated data at the intermediate 
selection level. 

At the final selection level, events with two or more isolated charged particles were 
required to have acollincarity and acoplanarity above 60°. The smaller of the two mo- 
menta had to be below 70% of the greater one, and below 10 GeV/c. For events with 
two reconstructed jets the scaled acoplanarity was required to be greater than 20°, and 
the acoplanarity and the acoUinearity greater than 60°. The acoplanarity distributions 
for the resulting samples of real and simulated data events are shown in figures 13(f) and 
14(f). 

5.5.5 Double tau cascade search 

The final state of X2X2 ~^ tttt differs from that of stau pair-production only in the 
presence of two soft taus from the stau decays. The search criteria for stau pairs (see 
section 5.2.3) could therefore be apphed with shght modifications, as follows. Instead of 
exactly two r-candidates, at least three were required, but the momentum of the third 
most energetic one had to be less than 4% of the beam momentum. The maximum 
number of tracks in the event was increased from seven to eight, and, as the W'''W~ 
background is negligible, the requirement on p^^^ versus Qm+ was removed. Since the 
only SM processes yielding final states with four taus arc ZZ, ZZ* or Z7*, the invariant 
mass of the two most energetic r-candidates was required to be below 40 GcV/c^, i.e. well 
below mi. Variables specific to the two-jet topology, e.g. acoplanarity, were evaluated 
using the two most energetic r candidates only, and the remaining requirements of the 
stau pair search were left unchanged. In the cases where these requirements differed for 
high and low AM, the less stringent selection was applied here. 

6 Results and limits 

In this section the results of the event selections and the estimated signal efficiencies 
are presented for each channel. No evidence of a significant excess with respect to the SM 
expectation was found in any of the channels, and limits on masses and cross-sections were 

set. In a relatively model-independent approach, cross-section upper limits are derived. 
Lower limits on the sparticle masses are also obtained, under assumptions which depend 
on the channel and will be specified case by case. All limits quoted are at 95% confidence 
level (CL). 

6.1 Limit computation 

Depending on the searches, limits were derived using the multichannel Bayesian 
method [38] or the modified frequentist hkelihood ratio method [33]. Both these ap- 
proaches allow search results in multiple channels to be combined, taking into account 
efficiency, expected background, number of candidates and centre-of-mass energy in each 
channel. 
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In the Bayesian multichannel approach [38] , the a posteriori probability density of the 
mean total number of signal events (assuming a uniform prior) is used to compute a 95% 
CL upper limit and compare it to the mean predicted for different signal assumptions. 
In the likelihood ratio approach, global multi-channel likelihood functions for the sig- 
nal plus background and background-only hypotheses are evaluated for the experimental 
outcome [33]. One then defines the confidence level in the signal plus background hy- 
pothesis as the probability, under this hypothesis, to obtain a lower likelihood ratio than 
experimentally observed. In the modified frequentist approach used here this confidence 
level is renormalised by the probability to obtain a value below the observed one in the 
absence of a signal. All points where the resulting confidence in the signal, CL^, is less 
than 5% are considered excluded. More details are given below channel by channel. 

The two methods typically give similar results. Both methods take a posteriori knowl- 
edge about the background into account, and thus give physically reasonable and con- 
servative results in the case of downward fluctuations of the background. Background 
fluctuations could nevertheless significantly affect the range of exclusion, and for this 
reason limits expected in the absence of a signal are also given below. 

6.2 Systematic uncertainties 

The excluded cross-section and mass ranges can be affected by the systematic uncer- 
tainties related to the SM background rate and signal efficiency determinations. These 
are mainly due to imperfections in: 

• the description of the differential cross-sections; 

• the modelling of fragmentation and hard gluon radiation; 

• the modelling of the detector response. 

The inaccuracy of the differential cross-section description for the various SM processes 
was studied extensively in the context of their measurement at LEP and is typically small, 
except for the hadronic 77 background. Modelling of fragmentation and hard gluon 
radiation was studied both in the context of background and efficiency. It mainly affects 
the searches for the stop quark, due to the uncertainties in the modelling of the stop 
fragmentation. However, the largest contribution to the systcmatics on both background 
rate and efficiency determination arises from the modelling of the detector response. 
This was addressed in several ways, depending on the type of the final state searched for. 
Below, a review of all of these systematic effects is given. The methods used to address 
them are then discussed. Finally, the methods used to propagate those uncertainties in 
the limit computation are briefly presented. 

6.2.1 Review of the main systematic sources 

Description of the differential cross-sections: 



The uncertainty on the differential cross-sections within the acceptance region of the 
most important SM processes is [25,39] less than 1% for W'''W~ production, around 
0.5% for Z7, 2% for ZZ and about 5% for Wez/g and Zc+e~. The uncertainty of 
the description of the Z7* process is typically 5%, reaching 20% in the region of very 
small masses of 7*. While the uncertainties are small for two- photon interactions giving 
leptons (less than 5%), the description of the hadronic 77 processes is subject to large 
uncertainties in some regions of the phase space. The effect was studied comparing various 

^"^This number is valid in the phase-space regions where these processes are not 77-Uke. 
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two-photon generators [39]. In the squark analysis (see section 5.3) at the preselection 
level PYTHIA [15] results were cross-checked with TWOGAM [21] and differences up to 15% 
were observed, with PYTHIA giving a higher background. It has been shown [40] that 
the PYTHIA generator represents a more correct background estimation. A systematic 
uncertainty of 15% on the hadronic 77 background was assumed. The impact of these 
uncertainties on the final background rate determination will depend on the importance 
of the two-photon background in the different cases. In most of the analyses presented 
in this paper the contribution of the hadronic 77 process at the final selection level is 
negligible. The search for multi-jet final states without isolated photons from heavier 
neutralino decays (see section 5.5) may serve as an example. In this case, 40% of the 
background originates from the W"'"W~ process, 40% from We Ve and ZZ processes, and 
15% from Z7. If the presence of photon(s) is required, 80% of the background originates 
from the Z7 process. In summary, the net uncertainties arising from the description of 
the above processes are of the order of a few percent at most. 

Fragmentation and hard gluon radiation modelling: 



The description of hard gluon radiation leads to uncertainties of at most 10% in the 
relevant parameter space, in particular for the rate of b-tagged four-jet events from the 
Z7 process. This is unimportant for most of analyses presented here. Fragmentation un- 
certainty related effects are typically of the order of 1%. However in the searches for the 
stop (see section 5.3), the stop hadronisation scheme is important and a dedicated genera- 
tor was used for these studies. The stop hadronisation was performed non-pcrturbatively 
and the e parameter of the Peterson function, which regulates the stop fragmentation and 
hadronisation, was varied. The analysis was applied to the different signal samples thus 
obtained and the results were compared. The relative systematic error on the efficiency 
as a consequence of imperfect signal simulation was taken to be 7% for AM> 10 GeV/c^ 
and 10% for AM< 10 GeV/c^. For the sbottom analysis, 7% was used for aU AM. 

Modelling of the detector response: 



The modelling of the detector response is the main source of systematics in the searches 
for new particles. The next section gives a review of the methods which were developed 
to address primarily the influence of this uncertainty on the SM background rate and 
efficiency determinations. 

6.2.2 Methods to evaluate detector response systematics 

The "re-weighting" method 



This method estimates the systematic error on the number of selected background events 
from imperfect simulation of the detector response, as well as from possible large un- 
certainties in the modelling of the differential cross-section for two-photon processes. It 
propagates the difference between the data and the simulated background at preselection 
level to the final selection level for the variables relevant for the analysis. 

The propagation is performed as follows. Each discriminating variable used in the 
event selection is histogrammed both for data and simulated background events at the 
preselection and final selection levels. From the preselection histograms a weight factor, 
-^DataZ-^Back, is Calculated for each bin containing at least 1 % of the total number of 
background events and a non-zero number of data events. For bins with zero data events. 
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it was checked that the SM expectation was in agreement with such an observation. 
These weight factors are then apphed, bin by bin, to the background histograms at the 
final selection level. The contribution to the systematic error of each variable is then 
determined by subtracting the total number of selected background events from the total 
number of re-weighted background events at the final selection level. 

The total positive (negative) systematic error of the number of background events is 
then computed as a quadratic sum over positive (negative) error contributions from the 
different variables: 



where N^^^{i) is the total number of re- weighted background events from variable i and 
-^Back is the number of selected background events without re-weighting. 

In this method the correlations between different variables are not taken into account, 
but, as the quadratic sum is performed over all error contributions, the total systematic 
error is overestimated. It was checked by the means of the Kolmogorov test that rescaling 
the distribution of any of the variables improves the agreement between the data and SM 
background for the other variables as well. Moreover, the statistical fiuctuation in the bins 
of the histograms used to calculate the weight factors lead to an additional ovcrcstimation 
of the error. Another drawback of this method is that the total systematic error obtained 
grows with the number of studied variables. 

The main advantages of this method are its simplicity and conservative error esti- 
mation. Because of its conservative nature, it was used in searches where the expected 
background is small, and background fiuctuations have bigger effect on the expected limit 
than the systematic error. For example, in the squark analysis the systematic error on 
the background estimated with this method was 15-60 % (see section 6.5.1). This had 
only a small effect on the final mass/cross-section limits. 

The "shaking" method 

The "shaking" method [41] attempts to correct at particle level the simulation description 
of the detector response and residual effects in the fragmentation. Changes in the particle 
multiplicity (by adding or removing particles) are made in the simulation, in such a way 
that the multiplicities of neutral and charged particles in bins of momentum and polar 
angle are well described for the hadronic data collected at the Z peak. The momentum 
and the polar angle of the added particles are randomly changed with respect to the 
parent particle, while staying in the same phase space region. 

The method was tested in the study of the systematics of the ZZ cross-section mea- 
surement [41]. The multiphcity adjustment is of the order of l%-2%. It leads to a better 
agreement between the data taken at energies above the Z peak and the SM background, 
typically predicting slightly higher observed missing energy and momentum. 

The method can be used to study the uncertainty on the background and efficiency in 
hadronic high multiplicity topologies, with or without isolated leptons (or photons). In 
these topologies, the method gives an estimation of the systematic uncertainty inchiding 
both the effect of "shaking" the variables and the consequences of additional reconstructed 
tracks which lie inside the cone for the isolated lepton or photon. 

This method was used to study the systematic error of background and efficiency in 
hadronic neutralino topologies, and the systematic error on the chargino (non-degenerate 
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scenario) detection efficiency in the hadronic, semileptonic and radiative topologies. Ta- 
bles 4 and 5 show, for chargino and neutralino searches respectively (sections 6.6.1 and 
6.7), the relative difference between the efficiencies obtained with and without "shaking". 
The values in table 4 correspond to high multiplicity hadronic and semi-leptonic final 
states and are averaged over several points with model parameters corresponding to the 
specific AM region (see section 5.4.1 for the AM region definition). The values shown in 
table 5 correspond to the different high multiplicity topologies in the neutralino search 
and they are averaged over the interesting range of masses. 



AM 


1 


2 


3 


4 


5 


6 


jets 


6.8 % 


8.8 % 


3.1 % 


0.5 % 


1.0 % 


3.4 % 




12.9 % 


4.1 % 


0.5 % 


0.4 % 


0.0 % 


0.9 % 



Table 4: The relative difference in the chargino detection efficiency obtained with and without "shaking" in 
the year 2000 simulation with the detector fully operational. The values shown correspond to the different AM 
regions of the hadronic (jets) and semi-leptonic topologies in the chargino production (see section 6.6.1). 



topology 


efficiency change 


multijets without 7 
multijets with 7 
cLCoplanar jets 


2.9 % 
1.3 % 
1.3 % 



Table 5: The relative difference in the neutralino detection efficiency obtained with and without "shafting" in 

the year 2000 simulation with the detector fully operational. The values shown correspond to the different high 
multiplicity topologies in the neutralino production (see section 6.7) and they are averaged over the interesting 
range of masses. 

Except for the lowest AM regions in the chargino analysis, where systematic uncer- 
tainties due to tracking and neutral energy reconstruction have the largest effect, the 
variation in the detection efficiency is very small. The efficiencies for the "shaken" signal 
are typically larger than for the "unshaken" signal. 

In hadronic neutralino topologies, the background estimated with the shaking method 
was typically 10% higher than the "unshaken" background. The variation of both back- 
ground and efficiency is consistent with "shaking" predicting slightly higher observed 
missing energy. 

Methods dedicated to low multiplicity topologies ("smearing" technique) 



The uncertainty on the efficiency of the muon, electron and photon identifications is 
expected to be a dominant effect in the low multiplicity topologies requiring identified 
particles. Studies on back-to-back di-muon events and back-to-back di-electron events 
and photons from radiative return to the Z peak point to background errors and rela- 
tive efficiency errors of the order of 3% each. The uncertainty on the track momentum 
reconstruction arising from the modelling of the detector response was also studied and 
a small effect was found, both on efficiency and background estimates. These estimates 
are relevant for neutralino and chargino leptonic topologies, and for all slepton searches. 



30 



6.2.3 Methods to propagate statistical and systematic uncertainties in the 
mass/cross-section limits 

Systematic and statistical uncertainties on the parameters involved in the calculation 
of the hmits (mainly, the uncertainties related to the SM background rate and signal 
efficiency determinations) were propagated into the final results using two methods. The 
analyses that rely on the modified likelihood ratio method adopted the procedure ex- 
plained in [33]. For the analyses that use the multichannel method of [38], a different 
procedure based on the same Bayesian approach was chosen [42]. A probability distribu- 
tion function (pdf ) was assumed for the efficiency and the background for every signal 
channel. Such a pdf was assumed to be cither gaussian or binomial, depending on the 
statistical accuracy of the estimations of the signal or background, and on the methods 
used to evaluate the uncertainty of the parameter. 

The effect of including systematics in the computation of the limits can become rel- 
evant whenever there are large (order of 10% or more) relative errors on the efficiency 
and, in particular, on the expected background. The degradation of the limits is much 
less significant, however, than would be indicated by the simplistic and over-conservative 
approach of reducing all efficiencies and increasing all backgrounds by one standard de- 
viation. 

6.3 LEPl limits 

In this section hmits on the masses of SUSY particles from LEPl data [43] are briefiy 
reviewed. In most cases the LEPl hmits have been superseded by LEP2 results, such 

as those presented in the present paper. However, for certain situations they arc still 
relevant. This is particularly true for limits deriving from comparisons of the measured 
Z decay widths to SM expectations. Such limits are relatively insensitive to the details 
of the decays of SUSY particles, although they depend on the couphng of the sparticles 
to the Z (which is affected by the sparticlc field composition). 

From 1990 to 1995 LEP was run at ccntre-of-mass energies near the Z resonance. 
Model independent fits to all the lineshape and asymmetry data have been carried out, 
giving accurate values of the resonance parameters [44]. The total decay width of the 
Z boson was measured with a precision of about 2.5 MeV/c^. Decay channels of the 
Z opened by new physics would increase the Z width. Thus the difference between the 
measured width and the SM value may be used to constrain SUSY models. If the new 
particles decay invisibly, limits can be derived in a straightforward manner from the 
comparison of the measured invisible width to the SM prediction. The combined LEP 
resuh gives Ff^^ < 2.0 MeV/c^ at 95% confidence level [45]. Whether the new particles 
are visible or invisible, they will contribute to the measured values of the total width F^. 
Confronting the measured Z width with the SM expectation an upper limit on the extra 
partial width F"^^ < 3.2 MeV/c^ was obtained [45] 

From the limit on PPf^, a limit on the sneutrino mass of 43.7 GeV/c^ may be ob- 
tained [45]. A lower mass limit for the lightest chargino of approximately 45 GeV/c^, 
independent of the field composition and of the decay modes, has been derived from the 
analysis of the Z width and decays. Limits for other sparticles depend both on masses and 
couphngs. To a large extent left-handed sleptons below 40 GeV/c^ can also be excluded 
using the agreement of the Z decay width with the SM prediction [45]. 

The composition of the neutralinos affects their production cross-sections and the 
light states may decouple from the Z. The production then proceeds through t-channel 

'^^In [46] a method designed to be more model-independent gives T""^" < 6.3 MeV/c^ using older data. 



31 



selectron exchange, giving a small cross-section if the selections are heavy. Hence no 
general limit on the LSP mass can be derived from LEPl. 

For the third family, the production rates may be affected by potential large mixing 
of the weak eigenstates. For the stau mixing angle giving the minimal cross-section, the 
coupling to the Z vanishes and no exclusion is possible using this method Also for 
squarks, the LEPl limits depend on the mixing angle. Left squarks below 45 GeV/c^ 
are excluded by the Z invisible width if they are nearly mass-degenerate with the LSP, 
decaying invisibly. If the decay is visible, the limit from the total width should be 
applied instead. Squarks with non-zero mixing cannot be excluded by this method, as 
the coupling to the Z varies and can vanish 

Direct searches at LEPl set mass hmits above 40 GeV/c^ for sfermions in the case of a 
decay into a fermion and a neutralino, with AM > 5 GeV/c^ and provided the production 
cross-section is not suprcssed due to t-channel contributions (in the case of selectrons) or 
to Z decoupling (in the case of third family sleptons and squarks) . 

6.4 Slepton searches 

6.4.1 Smuon and selectron seairches 

Efficiencies and selected events 

The efficiency for signal detection depends on the masses of the slepton and neutralino. 
The cuts used to reject the backgrounds resulted in typical signal efficiencies of 50% both 
for the selectron and the smuon channels. 

The number of events selected at each energy in the data, together with the estimate 
of the background, is shown in tables 6 and 7 for the selectron and smuon analyses. It 
can be seen that the principal background arises from leptonic decays of W pairs. 

The effect of systematic uncertainties on background and efficiency evaluation was 
studied with the "smearing" method (see section 6.2). The variations both in the detec- 
tion efficiency and in the background were found to be of the order of 3% on average. 





Selectrons 


Smuons 


year 


1998 


1999 


2000 


1998 


1999 


2000 


2-fermion events 


5.6 


8.4 


4.2 


2.1 


2.2 


1.4 


4-fermion events 


29.7 


40.9 


34.4 


18.8 


25.5 


20.8 


77 events 


1.7 


2.5 


2.2 


0.5 


0.6 


1.6 


Total 


37.0 


51.8 


40.8 


21.4 


28.3 


23.8 


Data 


40 


52 


49 


19 


23 


28 



Table 6l Selectron and smuon candidates in the different data sets, together with the number of background 
events expected. The systematic uncertainties on the background were estimated to be of the order of 3% (see 
section 6.2). 



Limits 

Selectron mass limits of the order of 50 GeV/c^ have been derived from single photon searches at early colliders, for a 
nearly massless photino LSP [47]. 

sufficiently light stop would contribute through loop corrections to the partial width F;,;,. 
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Sclcctron search 


(GeV) 


Preselection 


Intermediate 


Final 




Data 


MC 


Data 


MC 


Data 


MC 


189 


29460 


28558 ± 7 


1925 


1937.0 ±2.1 


40 


37.0 ± 0.4 


192 


4797 


4649 ± 7 


315 


319.7 ± 2.1 


6 


7.0 ± 0.4 


196 


13705 


14046 ± 9 


812 


847.2 ± 2.8 


21 


17.8 ± 0.6 


200 


14709 


14653 ± 8 


993 


1001.0 ± 2.5 


14 


18.5 ± 0.6 


202 


7125 


7154 ± 4 


474 


491.5 ± 1.4 


11 


8.6 ± 0.3 


205 


13229 


12687 ± 15 


832 


774.8 ± 4.7 


22 


14.7 ± 0.3 


207 


12735 


13315 ± 6 


810 


813.2 ± 2.0 


13 


15.4 ± 0.4 


208 


1202 


1231 ± 1 


78 


75.2 ± 0.2 


3 


1.4 ± 0.04 


206.5(*) 


8963 


9655 ± 6 


524 


495.6 ± 1.8 


11 


9.4 ± 0.5 
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205 
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8.5 ± 0.3 
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2854.7 ± 2.0 


11 


8.9 ± 0.4 


208 
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835 ± 1 
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0.8 ± 0.04 


206.5(*) 


6411 


6549 ± 6 


1696 


1621.6 ± 1.8 


10 


5.6 ± 0.5 


All 


70482 


70555 ± 2 


20626 


20707.5 ± 2 


70 


73.4 ± 0.4 



Table 7: Results of the selectron and smuon searches at the different selection levels and centre-of-mass energies. 
The number of events selected in data and expected from the SM simulation are given. Simulation errors are 
statistical; (*) indicates the 2000 data taken with the sector 6 of the TPC off. 

The results are presented in terms of excluded regions in the slepton-neutralino mass 
plane, obtained using all the analysed data sets combined with data taken previously at 
lower energies [2]. The method in [33] was used. 

Limits on slepton masses can be derived using several assumptions. In the MSSM, 
right-handed sleptons are expected to have lower masses and lower cross-sections for 
a given mass. Hence the assumption was made that only right-handed selectrons and 
smuons are sufficiently low in mass to be pair produced at LEP. This leads to conservative 
mass limits. 

Excluded regions in the slepton-neutralino mass plane were obtained taking into ac- 
count the signal efficiencies for each slepton-neutralino mass point, the cross-section for 
right-handed slepton production and the branching ratios squared for the direct decay 
£ — > £xij together with the number of data and background events kinematically com- 
patible with the mass combination under test. The estimate of the SUSY cross-section 
and branching ratios for each mass point were determined with the SUSY parameters 
tan P=1.5 and = — 200 GeV/c^. In addition, for the smuons a further exclusion curve 
was derived setting the branching ratio of /i ^ fiXi to 1. 

Figure 15 shows the excluded region for e^eR production. For a mass difference be- 
tween the selectron and the neutralino above 5 GeV/c^, right-handed selectrons are ex- 
cluded up to masses of 98 GeV/c^, for a neutralino mass up to 60 GeV/c^, beyond which 
the limit is weaker. For AM > 15 GeV/c^, the excluded mass range is up to 94 GeV/c^. 

Figure 16 shows the excluded regions for /xr/ir. The excluded region is shown both 
taking the branching ratios for each mass point with the SUSY parameters tan f3=1.5 
and /x = — 200GeV/c^ (lighter shaded region) and setting the branching ratio of /i — > i^Xi 
to 1 (darker shaded region). For the smuons the limit is determined at small neutralino 
masses, provided the mass difference between the smuon and the neutralino is above 
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5 GeV/c^. Assuming the branching ratio BR(/i —>■ /xx5)=l) masses up to 88 GeV/c^ are 
excluded. 

The effect of the systematic and statistical errors of background and efficiency on the 
mass limits was evaluated in a conservative way by changing the background and the 
efficiency by ± 3%. The resulting change in the limit was less than 1 GeV/c^. 

6.4.2 Stau search 

Efficiencies and selected events 

The efficiencies of the stau search have been determined using 5000 events for each 
point of a 1 GeV/c^ x 1 GeV/c^ grid in the (Mf,MLSp) plane, using SGV (see section 
4), and range from 20% to 30% for AM> 20 GeV/c^. The results have been verified with 
the full DELSIM detector simulation and analysis chain (figure 17(a)). In the low mass 
search, the efficiency only needed to be evaluated at a limited number of points, and in 
particular at AM= rur. Therefore, the full detector simulation could be used at AM= 
rur for Mf between 2 and 45 GeV/c^ (figure 17(b)). It was verified with SGV and, for a 
smaller number of Mf values, with DELSIM that the efficiencies were higher for higher 
AM. 

The systematic error on the efficiencies in the high AM region was obtained by re- 
versing the cuts designed to remove the r-pair and W^W" backgrounds, and comparing 
the number of selected events with the SM expectation. This method selected a sam- 
ple which contained 97.6% W'''W~ and r-pair events, while a possible signal would not 
exceed 1%. In terms of the most important kinematic distributions, the events selected 
were nevertheless quite similar to a high AM signal. The difference between data and the 
simulated SM processes was 5±5%, and 5% was taken as the estimate of the systematic 
uncertainty. For low AM, the uncertainty was estimated by the maximal scatter of the 
efficiency obtained with SGV with respect to the values obtained with DELSIM, and 
amounted to 15%. 

For the high mass stau-pair search, table 8 summarises the number of accepted events 
in the data together with the expected number of events from the different background 
channels. In all AM regions, good agreement with the SM expectation was observed. 
Most of the selected SM events in the simulation contained either one or two r's (17% 
and 67%, respectively). In two thirds of the events with less than two r's, the lepton 
mistakenly taken as coming from a r-decay had low momentum, i.e. was indistinguishable 
from a secondary lepton from a r-decay. In the low mass analysis, a total of 196 events 
kincmatically compatible with Mf = 25 GeV/c^ and AM = nir were selected, and the 
SM background was estimated to be 196. 1"*^^^°'^^. The contribution from 77 — rr to the 
background was 91%, and the remainder was other 77 processes (8%) and r-pairs (1%). 
In the very low mass analysis, 59 events compatible with Mf = 5.5 GeV/c^ and AM = 
m-r were selected, while the SM background was estimated to be 69.7 ± 1.8 events. The 
background was dominated by r-pairs (89%) and 77 processes (7%). The remaining 4% 
came from four-fermion processes. 

The systematic errors on the background estimates were obtained using the re- 
weighting method, as described in section 6.2.2. 

Limits 

Excluded regions in the (M^o,Mf) plane were derived, combining the analysed data 
with previous data sets [2]. For each Mf- M^sp mass combination, the predicted num- 



34 





Preselection 


Intermediate 


(GcV) 












Data 


MC 


Data 


MC 


189 


2949 


2916 ± 29 


80 


82 ± 2 


192 


473 


444 ± 4 


14 


14 ± 


196 


1265 


1333 ± 12 


36 


40 ± 1 


200 


1351 


1264 ± 11 


41 


39 ± 1 


202 


716 


652 ± 6 


12 


20 ± 1 


205 


1284 


1228 ± 11 


38 


34 ± 1 




1214 


1230 ± 11 


33 


36 ± 1 


208 


132 


116 ± 3 


3 


3 ± 


206.5(*) 


842 


857 ± 8 


28 


26 ± 1 


Total 


10226 


10040 ± 38 


285 


294 ± 3 


Channel 


Background composition 


w+w~ 




8 % 


85 % 


4-fcrmion 




% 


2 % 


2-fcrmion 




16 % 


1 % 


bhablia 




42 % 


% 


77 




34 % 


12 % 
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0.1 
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0.1 
0.1 
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0.1 



1.6 ± 0.5 
0.1 ± 0.0 
0.3 ± 0.0 
0.6 ± 0.1 
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0.4 ± 0.0 
0.5 ± 0.1 
0.2 ± 0.0 
0.5 ± 0.0 



FTTT- 

-2.3 



Total 



60 



58.4 ± 0.5 



36 



35.3 ± 0.5 



4.5 ± 0.5 



Channel 

w+w- 

4-fermion 
2-fermion 
bhabha 

77 



82 % 
3 % 
3 % 
% 
12 % 



Background composition 
77% 

4 % 

4 % 

% 
15 % 



54 % 
2 % 
% 
% 
44 % 



Table 8: Stau candidates, together with the number of background events expected, at the different selection 
levels and centre-of-mass energies. The column labelled "High AM" corresponds to the point with Mf=80 
GeV/c^ and Mlsp=0 GeV/c^ the one labelled "Medium AM" to Mf=80 GeV/c^ and Mlsp=40 GeV/c^ and 
the one labelled "Low AM" to Mf =65 GeV/c^ and Mlsp=60 GeV/c^. The composition of the SM background 
for all centre-of-mass energies summed is also given; (*) indicates the 2000 data taken with the sector 6 of the 
TPC off. 



ber of observed SUSY events was compared with the observed number of kinematically 
compatible events in data and simulated background. The method presented in [33] was 
used and systematic errors were taken into account when calculating CLs- The largest 
effect of the systematic uncertainty on the limit was observed at high AM, but never 
exceeded 800 MeV/c^. 

Figure 18 shows the tr excluded region. Figure 19 shows the excluded regions in the 
case of the mixing corresponding to the minimal production cross-section. The cross- 
section has a minimum at mixing angle 52° at LEP2 energies [48]. In addition, the 
limit for tl was also evaluated, and was found to be 84.7 GeV/c^ (expected hmit 84.9 
GeV/c^) for Mlsp — GeV/c^. The excess of candidates at low AM for high stau- mass 
seen in both figures is compatible with a statistical fluctuation: the observed limit was 
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everywhere contained inside the 95% confidence band for the expected hmit, calculated 
from the Poisson distributed background. 

Prom the low mass search, at the minimum cross-section, the lower limit on the stau 
mass was 26.3 GeV/c^ for AM = m^, 31.7 GeV/c^ for AM = 3 GeV/c^, and 40.5 GeV/c^ 
for AM = 4 GeV/c^. The corresponding expected limits were 26.3 GeV/c^, 35.9 GeV/c^, 
and 42.1 GeV/c^, respectively. The limit improves to 29.6 (31.1) GeV/c^ for tr (tl), 
with an expected hmit of 30.0 (31.9) GeV/c^, for AM — m-j-. The low mass search does 
not exclude below 6.3 GeV/c^ for AM below 3 GeV/c^. As shown in figure 20, to 
cover this region combination with the very low Mf search is needed. 

In summary, a stau mass limit can be set at 81.9 to 84.7 GeV/c^ (depending on 
mixing) for mass differences between the stau and the LSP above 15 GeV/c^. The same 
limits hold for LSP masses below 68 GeV/c^ and mass differences between the stau and 
the LSP above 6 GcV/c^. The expected limit in the same region ranges from 82.1 to 84.9 
GeV/c^. The lowest stau mass allowed is 26.3 GeV/c^ (any mixing-angle and any AM> 
nir)- 

6.5 Squark searches 

6.5.1 Efficiencies and selected events 

The efficiencies of the stop and sbottom signal selection are summarised in figure 21 
for the non-degenerate scenario. Systematic uncertainties on the efficiency determination 
have been estimated as explained in section 6.2.1. 

The number of events selected and the expected background at the final selection level 
are shown in table 9. The systematic uncertainties shown in the table have been estimated 
using the "re-weighting" method described in section 6.2.2. Moreover, in the mass limit 
computation, an additional relative uncertainty of 15% related to the determination of 
the hadronic 77 contribution (see section 6.2.1) has been added. 

In the nearly degenerate scenario, the efficiency for the stop signal is summarised in 
figure 22. The figure shows the variation of the efficiency for signal selection as a function 
of AM for different stop mass hypotheses, for centre-of-mass energies of 189 and 206 GeV. 
Some examples of the detection efficiency obtained with this analysis are quoted below: 
for a AM of 2 GeV/c^ and a stop mass of 70 GeV/c^, an efficiency of 2.4% is obtained at 
^/s = 189 GeV and of 4.8% at = 206 GeV. For a AM of 4 GeV/c^ and a stop mass 
of 80 GeV/c^, an efficiency of 4.0% is achieved at ^ = 189 GeV and of 7.7% at s/s ^ 
206 GeV. 

Table 10 gives the results for the nearly degenerate scenario in terms of number of 
events compared to the MC expectation after the preselection and at the final selection 
level. The systematic uncertainties were computed using the same method as for the non- 
degenerate case, and the same remark concerning the hadronic two-photon background 
applies. 

6.5.2 Limits 

Stop and sbottom cross-sections were calculated with the SUSYGEN program for 
two squark mixing angles. For purely left-handed squarks {9q — 0°), the cross-section 
is maximal. The squark mixing angle which corresponds to the Z decoupling is 56° 
for the stop and 68° for the sbottom, and it corresponds approximately to the minimal 
cross-section. The program ALRMC [33] was used and systematic errors were taken into 
account in the definition of CLg. 
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Table 9: Number of events selected by the squark analysis in the non-degenerate scenarios. 
The first errors are statistical and the second ones are systematic; (*) indicates the 2000 
data taken with the sector 6 of the TPC off. 

Figures 23 and 24 show the (Mq,M^o) regions excluded by the search for t — > cx5 and 

b — > bx5 decays, with the 100 % branching ratio assumption, both for purely left-handed 
states and for the states at the Z decoupling. Figure 25 shows the region excluded using 
only the nearly degenerate analysis for AM values between 2 and 20 GeV/c^. 

Table 11 shows the hmit on the squark masses as a function of AM obtained combining 
the two analyses. The introduction of the systematics in the confidence level calculation 
has no effect on these numbers. Stop masses lower than 71 GeV/c^ and sbottom masses 
lower than 76 GeV/c^ are excluded if AM > 2 GeV/c^ and AM > 7 GeV/c^ respectively, 
for any squark mixing angle. These hmits become 75 GeV/c^ and 93 GeV/c^ for purely 
left-handed squarks. 
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Stop (nearly degenerate) AM < 10 GeV/c^ 
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Table 10: Number of events selected in the stop analysis in the nearly degenerate scenario. 
The first errors are statistical and the second ones are systematic; (*) indicates 2000 data 
taken with sector 6 of the TPC off. 
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Table 11: Lower hmits on squark masses (in GeV/c^) as a function of AM from the 
squark analysis in the non-degenerate and nearly degenerate scenarios. 

6.6 Chargino searches 
6.6.1 Non-degenerate sceneirios 

Efficiencies and selected events 

The efficiencies of the chargino selection in the four topologies were computed sepa- 
rately for the 132 MSSM points using the C-ji function and the C-ji^^^^ of the corresponding 
topology and AM region. To pass from the efficiencies of the chargino selection in the 
four topologies to the efficiencies in the four decay channels, all the migration effects were 
computed for all the generated points of the signal simulation. Then the efficiencies of 
the selection in the four decay channels were interpolated in the (M-±,M^o) plane using 
the same method as in [41. These efficiencies as functions of M-± and M^o are shown in 
figure 26(27) for a mean centre-of-mass energy of 206(206.5) GeV with the TPC sector 6 
on(ofT). 
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Table 12 gives the total number of events selected in data and expected from the SM 
simulation after the preselections and after the final selections, for the different centre- 
of-mass energies. The total number of background events expected in the different mass 
windows and topologies is shown in tables 13, 14 and 15, together with the number of 
events selected. In all the topologies, the number of selected events in the real data is 
compatible with the expectation from the background simulation. 

The systematic errors shown in tables 13, 14 and 15 were obtained with the "re- 
weighting" method described in section 6.2.2. 
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Table 12: Numbers of events selected in data and expected from the SM simulation in the non-degenerate 
chargino search at the preselection and at the final selection level, for the different centre-of-mass energies collected 
during the years 1999 and 2000; (*) indicates the data collected with the TPC sector 6 off. Simulation errors are 
statistical. 
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Table 13: Number of events observed in data and expected number of background events in the different 
chargino search channels for all the data collected in 1999. The first errors are statistical and the second ones 
are systematic. The "re-weighting" method used to compute the systematics is described in section 6.2.1. 



To study the systematic effect on the detection efficiency both the "shaking" method 
and the "smearing" method were used (see section 6.2.2). 
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8 

7.7 ±0.9t2;| 


4 

5.6 ±0.7tli 


3 

2.5 ±0.5l J;J 


15 

16.3 ±1.3lJ;| 




25 < AM < 35 GeV/c^ 


Obs. events: 
Expect, events: 


1 

0.4 ±0.llg;f! 


13 

11.6 ±0.9+°i[ 


3 

4.0 ±0.3lJ;° 


3 

2.5 ±0.5t,V .1 


20 

18.5 ±l.ll2:i 




35 < AM < 50 GeV/c^ 


Obs. events: 
Expect, events: 


2 

2.4 ±OAto\ 


23 

26.8 il.Otj;! 


10 

8.0 ±0.4tg:° 


11 

10.3 ±0.5loJ 


46 

47.5 ±1.2l|| 




50 GeV/c^ < AM 


Obs. events: 
Expect, events: 


3.8 ±0.4t0-j 


40 

38.9 ±1.0ti:f 


22 

18.7 ±0.5to f, 


11 

10.3 ±0.5to g 


76 

71.7 ±1.3tf;^ 




TOTAL (logical .OR. between different AM windows) 


Obs. events: 
Expect, events: 


4 

5.0 ±0.5 


76 

69.7 ±2.2 


31 

29.0 ±1.3 


11 

10.9 ±0.7 


122 
114.6 ±2.7 



Table 14: Number of events observed in data and expected number of background events in the dif- 
ferent chargino search channels for all the events recorded in 2000 with the detector fully operational 
(< Ecm > = 206 GeV). The first errors are statistical and the second ones are systematic. The "re-weighting" 
method used to compute the systematics is described in section 6.2.1. 

For the high multipUcity topologies {jets & jji) the "shaking" method was used. The 
systematic uncertainty is larger in the two first AM regions, due to the more problematic 
event reconstruction. For these AM regions the relative systematic uncertainty on the 
detection efficiency determination is between 4% and 13%, but at larger AM the effect 
is less than 4%. The relative difference for these two topologies can be seen in table 4. 
All the results from the "shaking" method for these topologies gave a higher efficiency, 
indicating that the "unshaken" efficiencies can be regarded as conservative. 

For the leptonic {££) topology the "smearing" method was used. The results from 
this study did not give consistently higher efficiencies, but since the uncertainties are less 
than 1% for all the AM regions it showed that this topology is much less sensitive to 
systematic effects. 

In the case of the radiative topology, which consists of both high and low multiplicity 
events, both the "shaking" method and the "smearing" method were used to study the 
systematic uncertainty on the efficiency. The two methods gave compatible results in all 
the regions where they are both vahd and the uncertainty is less than 2%. 

Limits 

The simulated points were used to parameterise the efficiencies of the chargino selection 
criteria described in section 5.4.1 in terms of AM and the mass of the chargino. The 
values of AM, the chargino and neutralino masses, the cross-sections and the various 
decay branching ratios were then determined for a large number of points in the MSSM 
parameter space (//, M2, tan/3). From these values and the appropriate efficiencies, the 
number of expected signal events can be calculated. Taking into account the expected 
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Ecm = 206.5 GeV, C = 60.0 pb'^ 



Topology: 


jjt 


ee 


jets 


rad 


Total 




3 < AM < 5 GeV/c^ 


Obs. events: 
Expect, events: 




0.2 ±0.li;!;f, 


10 

18.8 ±1.3t!^;g 


2.1 ±0.4ir,.f, 


2 

1.3 ±0.1+1,1] 


14 

22.5 ±1.3l4;o 




5 < AM < 10 GcV/c^ 


Obs. events: 
Expect, events: 




0.2 ±0.lt°\l 


4 

3.6 ±0.5to° 


4 

3.9 ±0.7tlii 


2 

1.3 ±0.1+1-° 


10 

9.1 ±0.9lo;g 




10 < AM < 25 GeV/c^ 


Obs. events: 
Expect, events: 




0.2 ±0.ll°;J 


3 

2.31 ±0.3lo4 


7 

5.4 ±0.81^;^ 


2 

1.3 ±0.ltl-° 


12 

9.3 ±0.9lo 1 




25 < AM < 35 GeV/c^ 


Obs. events: 
Expect, events: 




0.4 ±0.ltoi 


3 

5.3 ±0.3t'^;!^ 


4 

3.2 ±0.3tJ;o 


2 

1.3 ±0.lli;» 


9 

10.3 ±0.5t2.2 




35 < AM < 50 GcV/c^ 


Obs. events: 
Expect, events: 




0.7 iO.llg ? 


10 

11.7 ±0.4l° ° 


6 

5.0 ±0.3tl l 


5 

4.3 ±0.2t J;3 


21 

21.6 ±0.6l5;2 




50 GeV/c^ < AM 


Obs. events: 
Expect, events: 


1 

2.0 ±0.llg;g 


16 

14.2 ±0.4t^;} 


9 

6.3 ±0.3to o 


5 

4.3 ±0.2lJ;o 


31 

26.8 ±0.5lt;| 




TOTAL (logical .OR. between different AM windows) 


Obs. events: 
Expect, events: 


1 

2.3 ±0.1 


31 

38.8 ±1.4 


16 

13.1 ±1.0 


5 

4.5 ±0.2 


53 

58.7 ±1.8 



Table 15: Number of events observed in data and expected number of background events in the different 
chargino search channels for all the events recorded in 2000 with the TPC sector 6 off (< Ecm > = 206.5 GeV). 
The first errors are statistical and the second ones are systematic. The "re-weighting" method used to compute 
the systematics is described in section 6.2.1. 

background and the number of observed events, the corresponding point in the MSSM 
parameter space can be excluded, if the number of expected signal events is greater 
than the upper limit on the number of observed events of the corresponding AM region, 
computed using the multichannel Bayesian formula (see section 6.1). Systematic errors 
were taken into account as described in section 6.2.3. 

Figure 28 shows the chargino production cross-sections as obtained in the MSSM at 
^/s = 208 GcV for different chargino masses for the large AM (AM > 10 GeV/c^) and 
low AM (AM = 3 GeV/c^) cases. The parameters M2 and ^ were varied randomly in 
the ranges GeV/c^ < M2 < 3000 GeV/c^ and -200 GeV/c^ < ^ < 200 GeV/c^ for 
three fixed different values of tan/?, namely 1, 1.5 and 35. The random generation of the 
parameters led to an accuracy on the mass hmit computation of the order of 50 MeV/c^. 
Two different cases were considered for the sneutrino mass: m,> > 1000 GeV/c^ (for 
AM > 10 GeV/c^) and m^>M^± (for AM = 3 GeV/c^). The mass limits are valid also 
for > 200 GeV/c^, since in this case both the efficiency of the chargino search and the 
chargino branching fractions are largely independent of if mi,> M-± 

The chargino mass limits are summarised in table 16. The dependence of the limit on 
AM and M2 assuming a heavy sneutrino is shown in figures 29 and 30. The behaviour 
of the curve in figure 29 depends very weakly on the relation between Mi and M2. Note 
that in figure 30, for a fixed large value of M2, the chargino mass limit is lower for positive 
// than for negative This is due to the higher degree of mass degeneracy between the 
lightest chargino and the LSP found for positive // compared to negative for a fixed 
value of M2. 

'^^This was investigated by a coarse scan of random parameter sets with |ju| > 200 GeV/c^, resulting in a mass limit 
greater than 103 GeV/c^. 
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For AM > 10 GeV/c^ with a large sneutrino mass (> 1000 GeV/c^), the lower limit 
on the chargino mass ranges between 102.7 GeV/c^ (for a mostly higgsino-like chargino) 
and 103.4 GeV/c^ (for a mostly wino-hke chargino). 

For AM = 3 GeV/c^ , the cross-section does not depend significantly on the sneutrino 
mass, since the region allowing small AM values is located where the chargino is higgsino- 
like, due to the assumption of gaugino mass unification. The lower limit for the chargino 
mass, shown in figure 28, is 97.1 GeV/c^. 



Case 


mo 




mtn 


N95% 




(GeV/c2) 


(GeV/c2) 


(Pb) 




< Ecm > = 192-208 GeV 


AM > 10 GcV/c2 


> 1000 


102.7 


0.39 


13.8 


AM = 5 GeV/c^ 


> 1000 


101.7 


0.57 


7.3 


AM = 3 GeV/c2 




97.1 


1.17 


18.1 



Table 16: Lower limits for the chargino mass, minimal pair-production cross-sections at 208 GeV for chargino 
masses below the limit and upper limits on the number of observed events. 

The results can be translated into a limit on the mass of the lightest neutralino [4] 
also shown in the (M~±,M^o) plane in figure 31. A lower limit of 38.2 GeV/c^ on the 
lightest neutralino mass is obtained, valid for tan /^ > 1 and a heavy sneutrino. This limit 
is reached for tan /3 = 1, = - 65.1 GeV/c^, M2 = 65.0 GeV/c^. 

6.6.2 Nearly mass-degenerate scenarios 

Efficiencies and selected events 

In the search for heavy stable charged particles, the three windows described in sec- 
tion 5.4.2 were searched for mass-degenerate charginos. No events were found in the 1999 
and 2000 data, where 0.51 ± 0.08 and 0.15 ± 0.03 events were expected, respectively. 
The efficiency for selecting a single chargino track by using this technique is shown in 
figure 32, as a function of the mass of the chargino and at the various centre-of-mass 
energies. 

In the search for kinks, 42000 chargino events with masses between 60 and 100 GeV/c^ 
and mean decay length of 50 cm were generated at the centre-of-mass energies between 
192 and 206 GeV. The events were used to map the selection and trigger efficiency for 
the single chargino, as a function of its mass and decay position. As an example, the 
efficiencies for a 75 GeV/c^ chargino at ^/s=196 GeV are plotted in figure 33. In the 
data of 1999 (2000 with full TPC, 2000 without the sector 6 of TPC), 5 (3, 1) events 
passed the selection, while 3.7± 1.0 (1.2 ±0.6, 0.5 ±0.2) were expected from the standard 
sources. 

In the search with the ISR photon tag a total of almost 3 million xtXi events was gen- 
erated at the centre-of-mass energies between 192 and 206 GeV; about 100 000 of them 
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had a high pi photon within the detector acceptance and were passed through the full 
detector simulation. The mass of the generated charginos ranged from 60 to 95 GeV/c^, 
and AM from 150 MeV/c^ to 3 GeV/c^. The selection efficiency was computed applying 
the selection cuts to the samples of simulated events. The trigger efficiency was param- 
eterised on the real data, separately for the isolated photon and for the system of few 
low momentum particles. Compton and Bhabha events were used to assess the trigger 
efficiency for single photons, as a function of the photon energy and polar angle. To 
estimate the trigger efficiency in events with few soft charged particles, parameterised in 
terms of the transverse momentum of the single particles and the total transverse energy, 
the redundancy of the trigger in several classes of two-photon events was used instead. 
The overall trigger efficiency for the whole event was finally considered as the convolution 
of the trigger efficiency of the single photon and that of the system of low momentum 
charged particles. The parameterisation was then applied on the simulation (both sig- 
nal and background events). Detection and trigger efficiency vary widely with Af-± and 
AM; they also depend on the field composition of the chargino, since the spectrum of 
the ISR radiation is different for higgsinos and gauginos. The efficiency (including the 
trigger efficiency) for a signal (higgsino) at 206 GeV is shown in figure 34. The small 
probability of radiating a photon with transverse energy Ej, above the chosen threshold 
is the main reason for having such a low efficiency. As the selection cuts varied across 
the plane {M^±, AM), so did the number of candidate events remaining in the data and 
the amount of background expected from the SM simulations available. Figure 35 shows, 
as an example, the number of events remaining in data and simulation for the year 2000 
sample with fully operational TPC, as a function of the position in the plane (M^±, AM). 

In this sample, when AM is below 1 GcV/c^ additional candidates are selected in the 
data but none in the simulated background. The number of events remaining after the 
dedicated preselection and at the final step, after the logical OR of all selections, together 
with the integrated luminosity used for the analysis, are given in table 17. The excess of 
events seen in the data is of the same kind as the excess observed in the analyses done at 
lower energies [6] and, as in the past, it can be qualitatively explained by the incomplete 
phase space coverage of the standard simulation used, and possibly some noise events 
like beam-gas collisions, as already discussed at the end of section 5.4.2. The excess at 
the preselection level has been verified to be almost entirely due to the cuts done at 
the generator level in the (hadronic) two-photon simulation used (in 1999 tighter cuts 
in the generation were applied to the hadronic 77 samples). After the final selection, 
the distribution of the 17 candidates in the data is compatible with the integrated lumi- 
nosities shown in table 17. Fewer candidates arc observed in the SM simulation than in 
the data. This is taken to indicate missing contributions in the modelling of two-photon 
background, or beam-gas interactions. 
Limits 

Having no evidence for a signal in any of the three searches for charginos nearly mass- 
dcgcncratc with the lightest ncutralino, regions in the plane (M^±, AM) can be excluded. 

First, the two searches for long-lived charginos were combined, assuming that in a XiXi 
event any of the two long-lived charginos can be tagged either as a kink or as a stable 
particle. Then the search with the ISR photon was considered for all events in which 
the chargino decay length was shorter than that required by the two other methods. In 
all cases, the data were combined with all previous data from the high energy phase of 
LEP [6]. 
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V3 
(GeV) 




Preselection 


Final scl. 


(pb-i) 


Data 


MC 


Data 


MC 


192 


25.4 


75 


40 ±2 





0.6± 0.1 


196 


76.2 


162 


118 ±4 


4 


1.3 ± 0.3 


200 


84.0 


201 


123 ±4 





1.4± 0.3 


202 


40.4 


109 


55 ±2 





0.6 ±0.1 


206 


163.0 


423 


393 ± 18 


11 


3.5 ± 1.3 


206.5(*) 


58.5 


130 


113 ±4 


2 


1.0 ± 0.4 


All 




1100 


842 ± 19 


17 


8.4 ± 1.5 



Table 17: Number of events selected in data and expected from the SM simulation in the ISR photon search 
for charginos nearly mass-degenerate with the Ughtest neutraUno. Numbers are given at the preselection and at 
the final selection level, for the different centre-of-mass energies. All year 2000 data were considered at the mean 
centre-of-mass energy of 206 GeV; (*) indicates the data collected with the TPC sector 6 off. Simulation errors 
are statistical. 

The limits obtained in this way are certainly model dependent: cross-sections, decay 
modes, and the spectrum of the ISR radiation itself all depend on the gauge composition 
of the chargino. They were obtained in the two SUSY scenarios in which a near mass- 
degeneracy between the lightest chargino and neutralino is possible: the M2 ^ /i scenario, 
in which the lightest chargino and neutralino are both almost pure higgsinos; and the 
/i ^ M2 scenario, in which the lightest chargino and neutralino are both almost pure 
gauginos. In the gaugino scenario, the gaugino mass unification at large scale must be 
violated in order to obtain low AM; in the higgsino scenario this is not mandatory, 
therefore this is the scenario to be taken into account in the constrained SUSY models. 
These limits arc shown in figure 36, separately for the different techniques. In the same 
figure, they are compared with the excluded region obtained in the search for larger AM 
charginos. 

With these new results, the AM-independent lower limit on the mass of the chargino 
becomes M^± > 75 GeV/c^ in the higgsino scenario and M-± > 70 GeV/c^ in the gaugino 
scenario with heavy sneutrinos. Both limits take into account a variation of tan f3 between 
1 and 50, and a variation of Mi, M2 and such that the mass difference between the 
chargino and the neutralino remains below 3 GeV/c^ and M2 < 2Mi < IOM2. In the 
higgsino scenario all sfermions are required to be heavier than the chargino, while in the 
gaugino scenario they must be heavy enough not to modify significantly the cross-section 
(only the sneutrino was required to be heavier than 500 GeV/c^) or the decay modes and 
widths. 

Uncertainties in the selection efficiencies and in the expected background were in- 
cluded in the limits in a Bayesian way as described in section 6.2.3, using the method 
in [42] . The most important contributions to the systematic error is the uncertainty in the 
background content: depending on the point, the excluded regions may vary as much as 
2-3 GeV/c^. Such an effect comes only from the uncertainty in the simulation available, 
since conservatively there was no attempt to increase the estimated background to take 
into account the regions of phase space not included in the simulation. Other systemat- 
ics, like the uncertainties in the trigger efficiencies, were similarly taken into account by 
considering the configuration which leads to the weakest limits. 
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6.7 Neutralino searches 

6.7.1 Efficiencies and selected events 

More than 1.2 million X1X2 events were simulated for different combinations of masses 
with M^o and Af^o ranging from 5 GeV/c^ to 100 GeV/c^ and from 10 GeV/c^ to 

200 GeV/c^, respectively, and for different X2 decay modes {qqXi-! /^^/^ Xi) e^e~Xi5 
fr ). The efficiencies were computed for each mass combination and parameterised as 
functions of the two masses. In addition, around 300 000 X2X3,4 events with cascade decays 
were simulated. For cascade processes, the efficiencies were parameterised as functions of 
M^o and a second parameter, chosen to be either M^o^—M^o when considering the decay 

modes X2 ~^ Xil X34 ^ > '^^ ^x%~^^X2 ^^^^ considering the decay modes 

X2 ~^ qqXi ^-iid X3 4 — > qqX2- The dependence on remaining parameters was found to 
be small in the relevant mass ranges and the efficiencies were averaged over them. 

The efficiency for X2X2 ~^ tttt was evaluated in a number of points with M^o between 
30 and 50 GeV/c^ and mass differences between the stau and the Xi ranging from 2 to 
5 GeV/c^. The X2 mass exceeded the Xi mass by 30 to 45 GeV/c^. These points were 
chosen because they correspond to a region in the CMSSM parameter space which cannot 
be covered by the chargino search nor, due to the small mass difference, by the stau search. 
Efficiencies between 10 and 20% were found, varying little with the chosen masses. These 
values were obtained with SGV and verified with DELSIM in a sub-set of points. 

The number of events selected at different centre-of-mass energies and selection levels, 
together with the expected SM background, are given in tables 18 to 22 for the different 
neutrahno search topologies discussed in section 5.5. The results of the sequential and 
likelihood ratio analyses for the acoplanar leptons and acoplanar jets were found to be 
comparable. The sequential analysis performed less well in the acoplanar jet search for 
large AM values, whereas their results were very similar in the acoplanar lepton channels 
and in the low AM region in general. 

Tables 22 and 23 summarise the results at the final selection level of all neutralino 
searches. Table 24 shows the main background sources contributing in each channel and 
the typical efficiency of each search for MSSM points where it is relevant. 

The effect of systematic uncertainties on background and efficiency evaluation was 
studied with the "shaking" method for high multiphcity topologies and with the "smear- 
ing" method for leptonic topologies (see section 6.2). The variations in the detection 
efficiencies were found to be small (the relative change was below 3% on average). The 
efficiencies for the "shaken" signal are typically larger than for the "unshaken" signal. 
The background estimated with the shaking method was typically 10% higher than the 
"unshaken" background. The variation of both background and efficiency is mostly due 
to a slightly higher observed missing energy in the "shaken" events. 



6.7.2 Limits 

In the absence of a signal, upper limits on neutralino production cross-sections were 
derived, using the dependence of the calculated efficiencies on the masses involved. The 
results obtained in different search topologies and at different centre-of-mass energies 
were combined using the Bayesian method [38] (see section 6.1). 

The limits for X1X2 production, as obtained from the searches for acoplanar leptons 
and jets, are shown in figure 37 assuming different branching ratios. 
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Acoplanar electrons search 


(GeV) 


Preselection 


Final 




Data 


MC 


Data 


MC 


192 


65 


60 ± 1 


5 


6.4 ± 0.1 


196 


135 


143 ± 3 


19 


15.7 ± 0.7 


200 


192 


180 ± 2 


12 


16.9 ± 1.0 


202 


68 


82 d= 1 


8 


9.7 ± 0.4 


205 


64 


77 ± 1 


18 


14.9 ± 0.7 


207 


72 


78 ± 1 


9 


15.0 ± 0.7 


208 


10 


6.5 ± 0.1 


4 


1.4 ± 0.1 


206.5(*) 


56 


62 ± 1 


16 


12.6 ± 0.7 




Acoplanar muons search 


(GcV) 


Preselection 


Final 




Data 


MC 


Data 


MC 


192 


102 


115 ± 3 


13 


7.7 ± 0.1 


196 


298 


289 ± 4 


18 


19.5 ± 0.4 


200 


340 


357 ± 3 


15 


21.0 ± 0.6 


202 


200 


160 ± 3 


14 


10.2 ± 0.2 


205 


179 


191 ± 2 


18 


19.5 ± 1.0 


207 


180 


193 ± 2 


20 


19.8 ± 1.1 


208 


24 


16 ± 0.4 


3 


1.9 ± 0.2 


206.5(*) 


179 


191 ± 2 


14 


16.3 ± 0.3 


V~s 




Acoplanar 


jets search 


(GeV) 


Preselection 


Final 




Data 


MC 


Data 


MC 


192 


927 


896 ± 6 


3 


3.1 ± 0.1 


196 


2191 


2218 ± 10 


13 


7.9 ± 0.3 


200 


2886 


2835 ± 12 


9 


10.4 ± 0.3 


202 


1263 


1250 ± 7 


7 


5.1 ± 0.2 


205 


1458 


1404 ± 8 


14 


12.7 ± 1.4 


207 


1451 


1416 ± 8 


15 


13.0 ± 1.4 


208 


133 


118 ± 2 


2 


1.2 ± 0.2 


206.5(*) 


1066 


1014 ± 8 


14 


7.8 ± 0.4 



Table 18: Results of the likelihood ratio acoplanar jets and acoplanar leptons searches at the different selection 
levels and centre-of-mass energies. The number of events selected in data and expected from the SM simulation 

are given. Simulation errors are statistical. The systematic uncertainties on the background were estimated to be 
of the order of 3% for the acoplanar leptons topologies and of the order of 10% for the acoplanar jets topologies 
on the final level (see section 6.2); (*) indicates the 2000 data taken with the sector 6 of the TPC off. 
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2.4 ± 0.3 
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1.2 ± 0.2 
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196 
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3.9 ± 0.5 
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9.3 ± 0.9 





4.3 ± 0.6 
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11.6 ± 0.4 
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4.3 ± 0.4 
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2.1 ± 0.3 
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5.5 ± 0.2 


205 
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6.3 ± 0.3 
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3.4 ± 0.2 
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11.5 ± 2.8 


207 
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7.4 ± 1.5 
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4.6 ± 1.2 


12 


10.9 ± 0.4 


208 


2 


0.7 ± 0.1 


1 


0.3 ± 0.1 


4 


0.9 ± 0.04 


206.5(*) 


6 


4.4 ± 0.6 


3 


2.5 ± 0.4 


10 


8.0 ± 0.5 


All 


41 


42.3 ± 2.1 


20 


22.3 ± 1.5 


53 


61.8 ± 3.0 



Table 19: Results of the acoplanar lepton and acoplanar jet sequential searches at the final selection level for 
the different flavours and centre-of-mass energies. The number of events selected in data and expected from the 
SM simulation are given. Simulation errors are statistical. The systematic uncertainties on the background were 
estimated to be of the order of 3% for the acoplanar leptons topologies and of the order of 10% for the acoplanar 
jets topologies (see section 6.2); (*) indicates the 2000 data taken with the sector 6 of the TPC off. 

Similarly, figures 38(a,b) show cross-section limits for X2Xi production {i= 3 or 4). 
Figure 38(a) shows the hmit obtained combining the results of the multijet and acoplanar 
jet searches in the case where X°— >X2Qq and X2~^X?qq. Figure 38(b) gives the corre- 
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Table 20: Results of the multijets without 7 and multijets with 7 searches at the diflFerent selection levels and 
centre-of-mass energies. The numbers of events selected in data and expected from the SM simulation are given. 
Simulation errors are statistical. The systematic uncertainties on the background were estimated to be of the 
order of 10% on the final level (see section 6.2); (*) indicates the 2000 data taken with the sector 6 of the TPC 
off. 
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Table 2 1 : Results of the multilepton and asymmetric tau searches at the different selection levels and cen- 
tre-of-mass energies. The numbers of events selected in data and expected from the SM simulation are given. 
Simulation errors are statistical. The systematic uncertainties on the background were estimated to be of the 
order of 3% on the final level (see section 6.2); (*) indicates the 2000 data taken with the sector 6 of the TPC 
off. 

sponding limits when x^^Xil^ obtained from the search for multijet events with a 
photon signature. 
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Table 22: Background and candidates in the X2X2 ~* tttt analysis. This analysis did not include data from 
the period in 2000 when one sector of the TPC was off. The total expected background over the three years 
was 1.10 ± 0.19 events, 55 % of which were two-fermion events, the remainder being four-fermion events. The 
systematic uncertainties on the background were estimated to be of the order of 5% on the final level (see section 
6.2.2). 
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Table 23: Results of the different neutralino seaxches. For any given search, events axe explicitly rejected if 
accepted by one of the searches appearing earlier in the table; (*) indicates the 2000 data taken with the sector 
6 of the TPC off. 

From the results of the search for X2X2 ~^ tttt described in section 5.2.3 , an equiv- 
alent production cross-section at ^/s = 206.7 GeV larger than 63 fb can be excluded. 
This hmit is valid for 30 GeV/(?< M^o < 50 GeV/c^ and Mf-M^o< 5 GeV/c^ in the 
gaugino region |//| >> M2, where this channel is important for constraining the LSP 
mass. To combine data at different centre-of-mass energies the cross-section dependence 
on ^/s given by CMSSM in the neutralino and stau mass ranges indicated above was 
assumed. The corresponding expected limit was 68 fb. 

The effect of the systematic and statistical errors of background and efficiency on the 
cross-section limits was evaluated using the method described in section 6.2.3 and found 
to be very small. 
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Table 24: The main backgrounds and the typical efficiency of all of the neutralino searches for MSSM points 
where it is relevant is shown. The efficiencies depend typically on the masses of the sparticles involved in the 
process. 

7 Combined exclusions and mass limits in the MSSM 
with gaugino and sfermion mass unification 

7.1 The method 

The method employed to set lower mass limits on the LSP mass and on the masses 
of other SUSY particles is to convert the negative results of the searches for charginos, 
neutralinos, sleptons and squarks into excluded regions in the {M2, fj.) plane for different 
tan (3 values, and then to find the minimal allowed sparticle masses as a function of tan (3. 
It is described in detail in [8,9]. The Higgs boson search in the maximal M^o scenario 
was used to exclude low values of tan f3 [11]. 

Unless stated otherwise, the limits presented in the following are valid for 
M2 < 1000 GeV/c^ and for the // region in which the lightest neutrahno is the LSP. 
The fi range depends on the values of tan/? and mo, and on the mixing parameters in 
the third family {At-, At, Ai,). Unless stated otherwise, for high values of itlq (above 500 
GeV/c^) the range between —2000 and 2000 GeV/c^ was scanned, but the scan range 
was increased if any limit point was found to be close to the scan boundary. 

7.1.1 Method of combining different secirches 

In the scan of the SUSY parameter space the efficiencies of the different searches, as 
obtained in the previous sections, were parameterised for the dominant channels, and 
used together with the information on the number of events selected in data and the 
expected number of background events. The excluded regions obtained with the different 
searches were then simply superimposed. 

In previous publications [8,9] it was verified that these parameterisations led to conser- 
vative results by comparing with a parallel approach, where these searches were applied 
to samples produced combining SGV with SUSYGEN (see section 4) to simulate simul- 
taneously all channels of chargino, neutralino, slepton and squark production and decay. 

The typical scan step size in and M2 was 1 GeV/c^ except in the region of the 
LSP mass limit, where the step size was decreased to 0.5 GeV/c^. The step size in mo 
was variable, the density of points being increased in regions of potentially difficult mass 
configurations. Special care was taken to set up the scan logics in such a way that no 
such configuration was overlooked. In particular, whenever two nearby scan points were 
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excluded by different searclies, tlie scan was performed witli smaller steps between these 
points to check the continuity of the exclusion. 

7.1.2 The influence of the mo value and of the mixing in the third family 

The unification of sfermion masses to a common mo at the GUT scale allows sfermion 
masses at the Electroweak Scale to be calculated as functions of tan/3, M2 and mo. In 
particular the masses of the sneutrino (u), the left-handed selectron and smuon {ei,,jli ), 
and the right-handed selectron and smuon (cr, /i^) can be expressed as 

(1) m£>2 ^ml + 0.77M| + 0.5mz^ cos 2/? 

(2) Ml^ml + 0.77M| - 0.27mz^ cos 2/3 

(3) M| = + 0.22M| - 0.23mz^ cos 2/3 

In the large mo scenario, mo = 1000 GeV/c^ was assumed, which implied sfermion 
masses of the same order. In this case only charginos, neutralinos and the Higgs bo- 
son could be produced at LEP and the limits arise from a combination of the searches 
for these particles described in the previous sections of this paper and in [11]. 

For large mo, the chargino pair-production cross-section is large and the chargino is 
excluded nearly up to the kinematic limit, provided M2 < 200 GeV/c^. 

It may also be remarked that at small M2, AM=M-± — M^o is large, resulting in 
increased background from W^W~ production. However, if is small as well, the 
chargino tends to decay via Xi xt^* the next-to- lightest neutralino which 
then decays by X2 ~^ Xi7 X2 ~^ Xi Z* ■ For setting the mass hmits, it is therefore 
important that the chargino search includes topologies with photons stemming from the 
decays xf ^ X2^* Xil^*i since the search for topologies with photons does not 
suffer from W^W^ background and is effective for large AM (close to mw)- 

Of the detectable neutrahno production channels {i.e. excluding XiXi): the X1X2 
XiX^ channels are dominant for large regions in the parameter space, but in order to cover 
as much as possible of the parameter space channels like X2X3 ^-iid X2X4 must also be 
considered, giving cascade decays with multiple jets or leptons in the final state. At large 
mo the production cross-section for all these neutralino production channels drops to very 
small values for |//| ^ 75 GeV/c^. This is because the two hghtest neutralinos have very 
small higgsino component (photino Xi and zino X2 ) their s-channel pair-production 
is therefore suppressed, while pair-production of heavier neutralinos is not kincmatically 
accessible. Nevertheless, for tan/3 < 1.5 and M2 > 68 GeV/c^ the neutralino exclusion 
reaches beyond the kinematic limit for chargino production at negative /i (see figure 39). 
This region of tan /3 is now also excluded by the searches for the production of the lightest 
Higgs boson [11]. 

For medium mo, 100 GeV/c^ ^ mo ^ 1000 GeV/c^, the X1X2 production cross- 
section in the gaugino- region (|//| 75 GeV/c^) grows quickly as mo falls, due to 
the rapidly rising contribution from the selectron i-channel exchange. Meanwhile, the 
chargino production cross-section in the gaugino region drops slowly, but it remains large 
enough to allow chargino exclusion nearly up to the kinematic limit for mo ^ 200 GeV/c^. 
For lower mo (~ 100 GeV/c^), the chargino production cross-section in the gaugino re- 
gion is close to its minimum, while the neutralino production cross-section is very much 
enhanced. Consequently, the region of the {fi, M2) parameter space excluded by the 

^^It is worth noting tiiat for tan/3 > 1 (tan/3 < 1) we have cos 2/3 < (cos 2/3 > 0), so tiic u is never heavier (Hghter) 
than the cl. 

^®If and/or mixing terms for the third family sfermions are sufficiently large, they can be light for large mo as well, 
this case will be discussed further on. 
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searches for neutralino production at small mo is larger than the one excluded by the 
search for chargino and neutralino production at large mo, as shown in [7-9]. 

For small mo, mo ^ 100 GeV/c^, and small M2, M2 ^ 200 GeV/c^, the situation is 
much more complicated because light sfermions affect not only the production cross- 
sections but also the decay patterns of charginos and neutralinos. Sleptons can also be 
searched for in direct pair-production. Excluded regions at small mo arise from the com- 
bination of searches for chargino, neutralino and slepton production. 

For small mo and M2 the sneutrino is light, and for M-± > m,> the chargino decay 
mode Xi ~^ is dominant, leading to an experimentally undetectable final state if M^± 
~ rriy. In the gaugino region, for every value of M2 and an mo can be found such that 
M-± ~ mp; therefore the search for charginos cannot be used to exclude regions in the 
(/i,M2) plane if very small mo values are allowed. The search for selectron production is 
used instead to put a limit on the sneutrino mass (and thus on the chargino mass), the 
selectron and the sneutrino masses being related by equations (l)-(3). The selectron pair- 
production cross-section is typically larger than the smuon pair-production cross-section, 
due to the contribution of ^-channel neutralino exchange. However, at |//| ^ 200 GeV/c^ 
the selectron production cross-section tends to be small and the exclusion arises mainly 
from the search for neutralino pair-production. 

Mixing between the left-handed and right-handed sfermion states can be important for 
the third family sfermions and can lead to light fi, bi and ti. Mass splitting terms at the 
Electroweak Scale proportional to mi-(^T— A*tan /?), mb{Ai,~ jitaji 13) , and mt{At—fi/ tan (3) 
were considered for f, b, and t respectively (see section 2.2). In the first instance 
Ar=At=Ab=0 was assumed, then the dependence of the results on Ar , Ai, and At was 
studied. These terms lead to fi, bi or ti being degenerate in mass with xi, or being the 
LSP for large values of The terms A-r,At, Ai, are arbitrary unless further constraints 
on the model are imposed. To derive a conservative limit on the LSP mass, which is valid 
for any stau mixing scenario, a model was used with no mixing in the sbottom or stop 
sector, but only in the stau sector. This maximises the region in the parameter space 
where the stau can be close in mass to the LSP. Moreover the trilinear coupling A^- was 
adjusted in every point of the parameter space to get the stau close in mass to the LSP 
mass. For the stau to be degenerate with the LSP in the interesting mass range, the 
corresponding mixing term has to be of the order of 12 000 GcV/c^. 

The results presented here are limited to the range of the /x parameter where the 
lightest neutralino is the LSP. 

7.2 Results 

7.2.1 The LSP mass limit for large mo 

From charginos searches alone a limit of 38.2 GeV/c^ on the lightest neutralino mass 
is obtained, valid for tan/3 > 1 and a heavy sneutrino {mo >300 GeV/c^). The hmit is 
reached for tan/3 = 1, = — 65.7 GeV/c^, M2 = 65.0 GeV/c^. This limit improves by 
~ 1 GeV/c^ due to the constraint from the search for neutralino production; thus from 
chargino and neutralino searches the LSP mass can be constrained to be M^o > 39.2 
GeV/c^, and the limit occurs at tan/3= 1. Figure 39 shows the region in the {11, M2) 
plane for tan 13—1 excluded by the chargino and neutrahno searches, relevant for the LSP 
mass limit at mo = 1000 GeV/c^. The lowest value of M^o not excluded by chargino and 
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neutralino searches occurs for tan/3=l, /j, = —76.5 GeV/c^ and M2 = 67.0 GeV/c^. For 
these parameters, X4X2 production and chargino pair-production are important. 

Figure 40 gives the lower hmit on M^o as a function of tan/3. At tan/3 ^ 1.2 the 
neutraUno search no longer contributes, the LSP limit is given by the chargino search, 
and its value reaches about half the limit on the chargino mass at large tan f3, where the 
isomass contours of xf and Xi iii the (/x,M2) plane become parallel. The rise of the LSP 
limit for small tan /3 can be explained by the change of the shape of these contours with 
tan/3, as illustrated in [8,9]. It should be noted that, because the chargino and neutralino 
masses are invariant under the exchange tan/3 ^ l/tan/3, the point tan/3 = 1 is the 
real minimum and the LSP limit for tan /3 < 1 can be obtained by replacing tan /3 with 
l/tan/3 in figure 40. 

For Ma < 1000 GeV/c^, Afii/tanP^Ve TeV/c^ (maximal M^o scenario used in [11]), 
and rrit — 174.3 GeV/c^, the tan/3 region 0.5< tan/3 < 2.36 is excluded by the Higgs 
searches. Thus including Higgs searches imposes a limit on M^o (see figure 40) of M^o > 
49.0 GeV/c^ for tan/3> 2.36 and M^o > 48.5 GeV/c^ for tan/3< 0.5. 

Thus the lightest neutrahno is constrained to have a mass: 

M^o > 49.0 GeV/c^ 

for mo = 1000 GeV/c^ and M2 < 1000 GcV/c^ and tan/3> 1.0. The limit occurs at 
the edge of the tan/3 exclusion given by the searches for the Higgs boson. However, if 
rrit = 179 GeV/c^, the tan/3 area excluded by searches for the Higgs boson shrinks to 
0.6< tan/3 < 2.0 and these limits worsen to M^o > 48.5 GeV/c^ for tan/3> 2.0 and M^o 
> 47.0 GeV/c^ for tan/3< 0.6. Thus the above limit worsens by 2 GeV if m* = 179 
GeV/c^ and the tan/3< 1.0 region is included. 

It should be noted that the tan /3 region excluded by Higgs boson searches is expected 
to shrink with the inclusion of complete one-loop corrections [49]. The LSP mass limit 
degrades in this case. Moreover, if there are sfermions lighter than Myfi/2, the hghtest 
Higgs boson will decay to them (see section 7.2.3). Such light sfermions are experimentally 
allowed only if they are degenerate in mass with the LSP. In this case, the lightest Higgs 
boson decays "invisibly". The mass limit M^o > 112.1 GeV/c^ set for an invisibly 
decaying Higgs boson [50] can be used to exclude small tan/3, and M^o > 48.5 GeV/c^ 
for tan /3> 2.0. 

7.2.2 The LSP mass limit for any mo 

Figure 40 also gives the lower limit on M^o as a function of tan /3 for any mo. The "any 
mo" limit resulting from chargino, neutralino and slepton searches follows the large mo 
limit up to tan /3= 1.4; then it increases more slowly due to the opening of the possibility 
of the chargino-sneutrino degeneracy, reaching 46 GeV/c^ at tan/3> 2.36 (edge of the 
region excluded by Higgs boson searches); finally, as discussed in more detail below, it 
falls to its lowest value, 45.5 GcV/c^, at tan/3> 5, due to small Mf-^-M^o, if mixing in 
the stau sector is of the form A^- — /xtan /3, and At- — 0. This limit is set by searches for 
XiXl and X2X2 production with xl fiX?- 

Thus 

M^o > 46.0 GeV/c2 

independent of mo, for tan/3 > 1 if there is no mixing in the third family [Ar = fitsmjS, 
Af, — //tan /3, At — /x/tan /3 ) , or if the mixing parameters leading to Mf^ —M^o < 6 GeV/c^ 
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are avoided. The Higgs boson search in the maximal M^o scenario was used to exclude 
0.5< tan/? < 2.36. 

LSP mass limits obtained with various assumptions are summarised in table 25. 

7.2.3 The dependence of the LSP limit on the mixing in the third family. 

Mixing in the third family affects in a complicated way the excluded regions obtained 
from Higgs, chargino, neutralino and squark searches, and a consistent discussion of the 
mixing is difficult. 

For example, the maximal M/^o scenario used to set tan/? limits from Higgs searches 
imphes A^ — ///tan /3 = ^/6 TeV/c^, thus a different At for every //. Such a scenario imphes 
as well that small mo values arc forbidden for M2 sufficiently small to be of interest from 
the point of view of the LSP mass limit Thus for consistency one should consider only 
the large mo scenario for the LSP limit (thus higher limit), whenever tan/? limits from 
the maximal M^o scenario are used. 

In the no-mixing scenario {Aj- = //tan/3. Ah — //tan/?. At = /i/tan/? ), the LSP 
limit occurs in the chargino-sneutrino degeneracy region, at tan/?> 1.5, where both the 
chargino and sneutrino mass limits are given by the selectron exclusion (see section 7.3). If 
there is no mixing in the stop sector, the tan /? region excluded by the Higgs searches grows 
to tan/?< 9.4, compared to 0.5< tan/? < 2.36 for the maximal M^o scenario. However, 
there is no improvement of the LSP limit as compared to Af^o > 46.0 GeV/c^ given 
above, as the tan/? dependence of the limit "flattens out" at large tan/? (sec figure 40). 

As discussed in previous papers [8,9], mixing in the stau sector can lead to a config- 
uration where Mf^ —M^o is small enough to make the fi undetectable and cause a blind 
spot both in fi and chargino exclusion. For At=AT-=Ai,=0, with the present data, the 
LSP mass limit occurs a) at large enough M2 that already for tan /?> 3 both M^^ and Mi^ 
are pushed above Mf^^^, and fi can become degenerate in mass with Xi^ and b) at large 
enough mo that selectron and sneutrino pair-production are not allowed by kinematics. 

For the above mixing assumptions, the mass of the lightest Higgs boson corresponding 
to the LSP limit point varies between 109 and 120 GcV/ c^, depending on the tan /? region. 
However, when the stau is light the lightest Higgs boson decays predominantly to fifi, 
thus "invisibly". The mass limit M^o > 112.1 GeV/c^ set for an invisibly decaying Higgs 
boson [50] can be used to exclude some of these points, but it is enough to change slightly 
the mixing in the stop sector (At) to push the Higgs mass above 112.1 GcV/cP. 

For such a mixing configuration and for tan/? ^ 3, the LSP limit is therefore set by 
the searches for stop and sbottom and those for X1X2 production with X2 ~^ ^1''") it 
falls to M^o > 45.5 GeV/c^ (see dot-dashed line in figure 40). This limit was verified to 
be robust when varying independently At^A^ and Af, in the range ± 1000 GeV/c^. 

As noted in previous papers [8,9], it is the appearance of the light sbottom and stop 
which protects the stau from being degenerate in mass with the LSP for very large mo 
values, where the X1X2 ^^'^ X2X2 production cross-sections are small. This is illustrated 
in figure 41, where the stop, sbottom, and stau masses at the largest allowed |/t| value at 

^"To avoid "tachyonic" mass solutions for the squarks and sleptons we must have mu + rrirr > (mu — rrirr)'^ + 4m^^ 
where m;^ is the off-diagonal mixing term, and mii,Tnrr arc the diagonal mass terms. For example, for the stop wc have: 
mir = mtop{At - /^/tan/3) and ma ~ mj; + 9M| + m?op + "'z 2/3(0.5 - 2/3sin^ 8w), rrirr ^ mf^ + 8.3M| + m^^^ + 
2/3rrt2 cos 2/3 sin^ ^Vl'; for an example value of At — /i/tan/3 = TcV/c^, the condition above sets a lower limit on a 
combination of and M|: mo^+S.SMg^ > 0.39 (TeV/c)^. Thus, if M2 < 190 GeV/c^ we must have mo > 300 GeV/c^. 

^^The "mixing -independent" (diagonal) terms of the mass matrices grow faster with M2 for squarks than for sleptons, 
and they have different dependence on tan/3. For example, for At=Ar=Ai,=0, /i=0, and tan/3=l, both the ti and bi are 
heavier than the fi; but they become lighter than the fi for large values. The mass hierarchy between fi, bi, and ti 
depends on M2 , tan (3, n, and mo • 
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the lowest non-excluded M2 are plotted. At large tan/5, the largest value occurs when 
the sbottom and stop masses have their lowest non-excluded (see section 6.5) values. In a 
pathological model where there is no mixing in the sbottom or stop sector {A}, — //tan /3, 
At — /i/tan/9) but only in the stau sector, fi can be made degenerate with M^o even at 
large values of mo and |/i| so that the x'j'x^ ^"^^ production cross-sections at LEP arc 
very small and the production of the Higgs boson and other sfermions is not accessible 
kinematically. However, for large ttiq values the XiXa cross-section is sufficiently large 
to set a limit on the production of "invisibly" decaying charginos (xf — > f\ u, fi— > xjr) 
from the search with an ISR photon (see section 6.6.2). The limit set on the neutralino 
mass in such a scenario is 39 GeV/c^ for tan/5>2 and is illustrated in figure 42. 

In mSUGRA (as defined in section 2), is in the range 3.3 m\i2—Q-^m?z < 1^^^ < 1^1 + 
3.8m1^2 tan/3>2 and a light stau cannot be degenerate with neutralino for large mo. 
Neutralino production cross-section is thus large, and neutrahno searches set a hmit on 
the LSP mass for small Mf^ — M^o which is close to the one obtained for heavy sleptons 
(about 50 GeV/c2). 



7.2.4 The influence of radiative corrections to the gaugino mass relations 
on the LSP limit. 

The relations between chargino, neutralino and gluino masses and and M2 are 
affected by radiative corrections of the order of 2%-20% [51]. However, only the rela- 
tive relations between chargino, neutralino and gluino masses are important from the 
experimental point of view, and here the corrections are much smaller. For example, the 
relation M-±{M-±)/MyO ~ 2 in the gaugino region, which is exploited to set a limit on 
the LSP mass, receives corrections only of the order of 2%. The ratio Mi/ai : M2/a2 can 
receive corrections of similar order [52], thus affecting M-±(M-±)/M^o ~ 2 in the same 
way. 

It is enough to lower the LSP limits presented here by 1 GeV/c^ (2%) to conservatively 
account for the 2% change of the M-±(M-±)/M^o ratio. 



7.3 x-^ mass limits for any mo 

Figure 43 shows the chargino mass limit as a function of tan/? for M2 < 200 GeV/c^. 
The lowest non-excluded chargino mass is found at MSSM points very close to those 
giving the LSP mass limit, and the arguments presented in section 7.2.2 also apply to 
explain the dependence of the chargino mass limit on tan/3. For tan/3 ^ 1.4 the limit 
occurs at large mo values. For 1.4 ^ tan/3 ^ 3 and M2 < 200 GeV/c^, the limit occurs 
at small mo in the chargino-sneutrino degeneracy region. It falls at tan /3 ^ 4 because 

of the small AM = - M^o. 

' Xi 

The lightest chargino is constrained to have a mass: 

M-± > 94 GeV/c^, 

independently of mo, for tan/3 < 40, M2 < 1000 GeV/c^, if either there is no mixing in 
the third family {A^- — //tan/3, Af, — /(tan/3, A^ — /(/tan/3) or the mixing parameters 
leading to Mf^ — M^o < 6 GeV/c^ are avoided. 

If mixing in the stau sector is of the form Aj. — /itan/3, and Aj. = 0, the limit falls to 
90 GeV/c^, at tan/3> 3, due to small Mf^—M^o. This limit is robust when varying At^A-j- 
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and Ah independently in the range ± 1000 GeV/c^ In the "arbitrary" mixing scenario 
described above, where there is no mixing in the stop and sbottom sector, but only in 
the stau sector, this hmit falls to 80 GeV/c^ and it is independent of the stau mixing. 

The chargino mass limits for large M2 values, where the chargino can be degenerate 
in mass with the LSP are close to 75 GeV/c^ (sec section 6.6.2). 

Chargino mass limits obtained with various assumptions are summarised in table 25. 

7.4 i> and gr mass limits for any mo 

The sneutrino and the cr mass limits are: 

> 94 GeV/c^ and M~^^ > 94 GeV/c^. 

These limits, shown in figure 43, were obtained assuming no mass splitting in the third 
sfermion family (Ar — //tan/3), implying Me^—Mj^—Mf^—Mfi^, as this gives the lowest 
values. 

These limits result from the combination of slepton and neutralino searches. The 
selectron mass limit (see figure 43, dotted curve) is vahd for —1000 GeV/c^ < A* < 1000 
GeV/c^ and 1 < tan/3 < 40 provided that M^^ - M^o > 10 GeV/c^ and it allows a 
limit to be set on the sneutrino mass as shown in figure 43 (dashed curve). The sneutrino 
mass limit is expected to rise for tan/3 < 1, the sneutrino being heavier than the e^ for 
small tan/3. The selectron mass limit for tan/3 = 1.5 and /i = —200 GeV/c^ is presented 
in section 6.4.1. 

Slepton mass limits obtained with different assumptions are summarised in table 25. 

7.5 MSSM parameter space exclusion 

The limits on Higgs, chargino, neutralino, slepton and squark production used in the 
previous sections to set a mass limit on the LSP and other particles can be also used to 
exclude regions in a parameter space of the CMSSM. The excluded regions in the (/i, M2) 
plane for tan /3 = 35 and two assumptions about mo and mixing in the third family are 
shown in figure 44. 

8 Summary 

Searches for charginos, neutralinos, sleptons and squarks in c+c~ collisions at ccntre- 
of-mass energies up to 208 GeV were performed with the DELPHI detector at LEP. 
No evidence for a signal was found in any of the channels and 95% CL upper limits 
on the production cross-sections were derived. Under assumptions that depend on the 
channel lower limits on the masses of SUSY particles were set. In particular, in the 
framework of constrained MSSM scenarios with gravity-induced SUSY breaking, regions 
of the parameter space can be excluded, and these exclusions can be translated into limits 
on the masses of SUSY particles. The combination of the results of the different search 
channels is crucial to ensure the best possible coverage of the parameter space. 

A summary of the obtained lower limits on the masses of SUSY particles is shown in 
table 25. The results presented extend and confirm previous exclusions set by DELPHI 
(see [2] to [9]) and by the other LEP experiments [53]. 
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Particle 


Validity conditions 


Mass limit 
(GcV/c2) 


eR 


tan/?=1.5, /x=-200, AM>15 


94 




CMSSM, AM>10 


94 




BR(/i^ AiX°) = l , AM>5 


88 




CMSSM, AM>10 


94 


f 


BR(f-^ rx°)=l, AM> 


26 




BR(f— rx'^)=l, AM>15, no mixing 


85 




BR(f^ Ty^)=l, AM>15, minimal cross-section 


82 




CMSSM, (Mej,-M^o)>10 


94 


b 


BR(b^ b x°)=l,AM>7, no mixing 


93 




BR(b— > b x°)=l, AM>7, minimal cross-section 


76 




BR(t^ c x°)=l, AM>10, no mixing 


96 


t 


BR(t^ c AM>2, no mixing 


75 




BR(t— > c x'^)=l, AM>10, minimal cross-section 


92 




BR(t— > c x'^)=l, AM>2, minimal cross-section 


71 




mf^>1000, AM>10, Ml =~ O.5M2, 


102.7 




M, > M AM>3 


97 




M, > ,1/ anv AM. ,1/^ =- 0.5,1/.; 


75 




/;/,>> 300. // > Mo. IK) gaugiiio mass unification, any AM 


70 




CMSSM. AM>3. any ni,-,. no nuxing or AM(f-v")>6 


94 




CMSSM. any nio. any Mo. tan i <10. mixing A^=A(;=A,=0 


90 




CMSSM, high mo, tan/3 >1, maximal mixing in t sector 


49 




CMSSM, any mo, tan/3 <40 no mixing or AM(f-x°)>6 


46 




CMSSM, any mo, tan/3 <40, mixing A^^A^^Aj = 


46 




CMSSM, any mo, Ktan/3 <40, mix. A.,=A6=0, A* = \/6 TeV/c^ 


49 



Table 25: Summary of mass limits for supersymmetric particles and their validity conditions. In each line of 
the table AM is the mass difference between the corresponding sparticle and the LSP. All masses and AM values 
are in GeV/c^. CMSSM refers to a model with gauge and sfermion mass unification, where /^ however is a free 
parameter (see section 2). Neutralino mass limits should be lowered by 1 Ge\f /<? if the radiative corrections 
of [51] are taken into account. 
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Figure 1: Comparison of data and simulation in the selectron channel at preselection 
level. The dots with error bars show the data, shaded histograms show the simula- 
tion. Plots include data taken in the year 2000 when the DELPHI detector was fully 
operational. The plots show: (a) the acoplanarity, (b) the transverse momentum, (c) 
the opening angle, (d) the momentum of the leading charged particle. Possible sig- 
nals corresponding to the mass combinations Me=90 GeV/c2,M^o=10 GeV/c^ (sohd) 
and Me=50 GeV/c^,M^o=40 GeV/c^ (dashed) are shown by the superimposed open his- 
tograms. The signal normalisation is arbitrary. 
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Figure 2: Comparison of data and simulation in the smuon channel at preselection 
level. The dots with error bars show the data, shaded histograms show the simula- 
tion. Plots include data taken in the year 2000 when the DELPHI detector was fully 
operational . The plots show: (a) the acoplanarity, (b) the transverse momentum, (c) 
the opening angle, (d) the momentum of the leading charged particle. Possible sig- 
nals corresponding to the mass combinations M^=90 GeV/c^,M^o=10 GeV/c^ (sohd) 
and M^=50 GeV/c^,M^o=40 GeV/c^ (dashed) are shown by the superimposed open his- 
tograms. The signal normalisation is arbitrary. 
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Figure 3: A preselection comparison of data and simulation in the stau analysis. The plots 
show: (a) the number of charged particles, (b) the square of the transverse momentum 
with respect to the thrust axis, (c) the acoplanarity, (d) the missing transverse momentum 
divided by the maximum missing transverse momentum in two-photon events with no 
beam-remnant electrons in the detector acceptance (ie. in "no-tag" events). The dots 
with error bars show the data, while the simulation is shown shaded. A typical signal 
{Mf= 83 GeV/c^, Mlsp = GeV/c^) is shown by the superimposed open histogram, 
with arbitrary normalisation. 
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Figure 4: Comparison of data and simulation at the preselection level in the 
non-degenerate squark analysis. Plots include all DELPHI data from 189 to 208 GeV: (a) 
visible energy, (b) acoplanarity, (c) combined b-tagging variable, (d) maximal transverse 
momentum of a charged particle. The expected signal distributions at a/s = 200 GeV 
are shown for one possible stop and sbottom signal (Mq=90 GeV/c^, M^o = 60 GeV/c^), 
with arbitrary normalisation. 
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Figure 5: Number of events as a function of the signal detection efficiencies in 
the non-degenerate squark analysis: (a) 1998 data at 189 GeV: stop analysis 
for AM > 20 GeV/c^, (b) 1999 data from 192 to 202 GeV: sbottom analysis 
for AM > 20 GeV/c^ (c) 2000 data with TPC sector 6 off: stop analysis for 
5 < AM < 20 GeV/c^, (d) 2000 data with TPC sector 6 on: sbottom analysis for 
5 < AM < 20 GeV/c^. The efficiencies are for a given combination of squark and LSP 
masses, indicated in the parentheses. 
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Figure 6: Comparison of data and simulation in the nearly degenerate squark analysis at 
preselection and final selection level. The distributions of (a) the total transverse energy, 
(b) total number of charged particles, (c) angle between the missing momentum and the 
beam axis and (d) transverse momentum are shown at preselection level. At the final 
selection level the distributions of (e) the acoplanarity and (f) the transverse momentum 
are shown. All DELPHI data from 189 to 208 GeV is included. A signal corresponding 
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Figure 7: a), b) and c) show comparisons between real data (dots) and simulated back- 
ground events (histogram) for the jets, ii and rad topologies respectively, using a logical 
OR of the 6 preselection cut functions of the corresponding topology. The dashed lines 
indicate how a characteristic chargino signal would appear (arbitrary normalisation), d) 
shows the number of events selected by the standard chargino analysis as a function of 
the Cti^j^^ cut in the jjl topology for 25 < AM < 35 GeV/c^. The squares are the data 
and the solid line is the background simulation. The dotted curve shows a possible signal, 
M-±=102.8 GeV/c^ M^o=73 GeV/c^ of 0.5 pb. In all cases the data collected in the 
year 2000 with the TPC sector 6 on are shown. 
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DELPHI non-degenerate xtxi search 
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Figure 8: a), b) and c) show comparisons between real data (dots) and simulated back- 
ground events (histogram) for the jets, ii and rad topologies respectively, using a logical 
OR of the 6 preselection cut functions of the corresponding topology. The dashed lines 
indicate how a characteristic chargino signal would appear (arbitrary normalisation), d) 
shows the number of events selected by the standard chargino analysis as a function of 
the Cti^j^^ cut in the jjl topology for 25 < AM < 35 GeV/c^. The squares are the data 
and the solid line is the background simulation. The dotted curve shows a possible signal, 
M-±=102.8 GeV/c^ M^o=73 GeV/c^ of 0.5 pb. In all cases the data collected in the 
year 2000 with the TPC sector 6 off are shown. 
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Figure 9: Some of the variables used in the selection for mass-degenerate charginos with 
an ISR photon tag. In the plots on the left the data (dots) are compared with the SM 
expectations. On the right, as an example, the corresponding distributions (with arbitrary 
normalisation) are shown for the signal with M{xt) = 80 GeV/c^ and AM = 1 GeV/c^. 
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Figure 10: On the left the comparison between the real and simulated data is shown at 
preselection level for different variables for the qq, ee, and fifi topologies in the neutralino 
search. On the right the visible energy distribution is shown for the events selected at the 
final stage, after the likelihood selection. The signal distribution contains the contribution 
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Figure 11: Comparison of real data from 1999 {^/s in the range 192-202 GeV) and 
simulation for the neutralino multijet selection without photons (a,c,e) and with photons 
(b,d,f), at the preselection (a,b), intermediate (c,d), and final selection(e,f) stages. The 
distributions expected for X3X2 production with Xg X^*!*! X2 decaying to x^qq 
or Xil are shown as dashed histograms where M^o = 112 GeV/c^, M^o=75 GeV/c^ and 
M^o=41 GeV/c^. The signals are normalised to cross-sections of 0.8 pb (e) and 0.1 pb 
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Figure 12: Comparison of real data from 2000 (y^ in the range 204-208 GeV) and 
simulation for the neutralino multijet selection without photons (a,c,e,g) and with photons 
(b,d,f,h), at the preselection (a,b), intermediate (c,d), and final selection(e,f,g,h) stages. 
In (g,h) only the data collected with the TPC sector 6 on are shown. The distributions 
expected for X3X2 production with X3 s-^d X2 decaying to x'Jqq or x'Jt are shown 

as dashed histograms where M^o=112 GeV/c^, M^o=75 GeV/c^ and M^o=41 GeV/c^. 
The signals are normahsed cross-sections of 0.8 pb (e), 0.1 pb (f,h) and 0.4 pb (g). 
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Figure 13: Comparison of real data from 1999 {^/s in the range 192-202 GeV) and simu- 
lation for the neutralino multilepton search (a,c,e) and asymmetric tau search (b,d,f), at 
the preselection (a,b), intermediate (c,d), and final selection (e,f) stages. The dashed line 
in (e) shows the multilepton signal expected for X3X2 production with Xs X2 and 
xl decaying to where M^o=103 GeV/c^, M^o=51 GeV/c^ and M^o=45 GeV/c^. 

In (f) the dashed hne shows the tau cascade signal for X2X1 production with X2 ~^ ^'^ ^-nd 
f ^ X?r where M^o=75 GeV/c^ M^=44 GeV/c^ and M^o=37.5 GeV/c^. The signals 
are normalised to a cross-section of 0.4 pb. 
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Figure 14: Comparison of real data from 2000 {^/s in the range 204-208 GeV) and simu- 
lation for the neutralino multilepton search (a,c,e) and asymmetric tau search (b,d,f), at 
the preselection (a,b), intermediate (c,d), and final selection (e,f) stages. The dashed line 
in (e) shows the multilepton signal expected for X3X2 production with Xs X2 and 
xl decaying to where M^o=103 GeV/c^, M^o=51 GeV/c^ and M^o=45 GeV/c^. 

In (f) the dashed hne shows the tau cascade signal for X2X1 production with X2 ~^ ^'^ ^-nd 
f ^ X?r where M^o=75 GeV/c^ M^=44 GeV/c^ and M^o=37.5 GeV/c^. The signals 
are normalised to a cross-section of 0.4 pb. 
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Figure 15: Excluded region in the (Mgj^ Mlsp) plane. The solid line shows the expected 
limit and the shaded region shows the obtained limit. The cross-section and the branching 
ratios for each mass point were determined with the SUSY parameters tan /3=1.5 and 
/i = — 200GeV/c^ (for a discussion on the low selectron mass region see section 6.3). 
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Figure 16: Excluded region in the (M^j^ Mlsp) plane. The shaded region shows the 
obtained limit. The lighter shaded region is excluded taking the branching ratios for 
each mass point with the SUSY parameters tan /3=1.5 and yU = — 200 GeV/c^. The solid 
line shows the corresponding expected limit. If the branching ratio of ft —>■ fiXi is set 
to 1, the darker shaded region is also excluded (for a discussion on the low smuon mass 
region see section 6.3). 
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Figure 17: The selection efficiency in the stau analysis. The plots show: (a) Efficiency as 
a function of Mlsp at Mf = 80 GeV/c^. The line shows the result of SGV, and the points 
with error bars show the result of the full DELSIM simulation. The shaded area indicates 
the statistical uncertainty of the estimate, and the vertical line shows the position of the 
transition between the two AM-regions. (b) Efficiency as a function of Mf at AM= m,- 
for the low mass analysis (filled circles) and the very low mass analysis (open circles). 
The points show the results of the full simulation The lines are polynomial ffis to the 
points. 
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Figure 18: Excluded region in the (Ms- , M^sp) plane. The shaded region shows the 
region excluded if r— > r LSP is the only allowed channel, and the thin solid line shows the 
corresponding expected limit. The dotted lines show the region excluded if the branching 
ratio of f— > r LSP has the values indicated in percent next to the lines, and it is assumed 
that the analysis has no sensitivity to other decay-modes. The dash-dotted line indicates 
AM = 171^-. The observed limit is everywhere within the 95% CL band for the expected 
limit shown as the hatched area. 
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Figure 19: Excluded region in the (M^, Mlsp) plane obtained for the minimal f 
pair-production cross-section. The shaded region shows the region excluded if r 
LSP is the only allowed channel, and the thin solid line shows the corresponding ex- 
pected limit. The dotted lines show the region excluded if the branching ratio of r 
LSP has the values indicated in percent next to the lines, and it is assumed that the 
analysis has no sensitivity to other decay-modes. The dash-dotted line indicates AM = 
rrir- The observed limit is everywhere within the 95% CL band for the expected limit 
shown as the hatched area. 
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Figure 20: The value of CLg versus Mf obtained for the minimal f pair-production 
cross-section, and AM = rrir, which corresponds to the weakest limit. The solid curve 
shows the obtained CLg, and the dashed line shows the corresponding expected value. 
The two sets of curves correspond to the "very low mass analysis" and the "low mass 
analysis" , as indicated in the figure. 
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Figure 21: Squark signal detection efficiencies (in %) in the plane (Mq^,M^o) in the 
non-degenerate scenario. 
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Figure 22: Signal efficiency of the stop search in the nearly degenerate scenario for dif- 
ferent mass hypothesis as a function of AM = M^-^ — M^o, for ^/s=189 and 206 GeV. 
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Figure 23: Excluded regions for the stop search in the plane (Mj^,M^o) for purely left 
handed stops (top) and for the states at the Z decoupling (bottom). The shaded areas 
show the observed excluded regions and the lines correspond to the expected exclusions 
(for a discussion on the low stop mass region see section 6.3). 
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Figure 24: Excluded regions at for the sbottom search in the plane (Mj^^,M^o) for purely 
left handed sbottoms (top) and for the states at the Z decoupling (bottom). The shaded 
areas show the observed excluded regions and the lines correspond to the expected ex- 
clusions (for a discussion on the low sbottom mass region see section 6.3). 
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DELPHI 189 to 208 GeV at 95% CL 
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Figure 25: Stop mass limit as a function of AM for a purely left-handed stop and for the 
stop mixing angle which minimises the cross-section. 
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Figure 26: Chargino pair-production detection efficiencies (%) for the four decay channels 
a) b) jefo, c) a and d) rad, at 206 GeV in the {M^±^M^o) plane. TPC sector 6 is 
on and a stable Xi is assumed. The shaded areas are disallowed in the MSSM scheme. 
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DELPHI xlxi efficiencies at 206.5 GeV 



e^e' — > ~> jjl detection efficiency 



e^e" — > ~> js'*^ detection efficiency 




80 100 

M(x^) (GeV/c') 

e^e — > X*X ~> U detection efficiency 



80 100 

M(x^) (GeV/c') 

e^e" — > ~> YY^ detection efficiency 



> 100 

o 

°§ 80 
60 
40 
20 




; Stable x? 




c) 








--^TO/ / >40 






1 \ >30 ^^^^ 






: >2o ) i 

■ 1 1 ' 1 » 


, ,> 1 , , , , 1 





60 80 100 

Mix*) (GeV/c') 




80 100 

M(x*) (GeV/c') 



Figure 27: Chargino pair-production detection efficiencies (%) for the four decay channels 
a) jj^, b) jets, c) and d) rac?, at 206.5 GeV in the (M-j^±,M^o) plane. TPC sector 6 is 
off and a stable Xi is assumed. The shaded disallowed in the MSSM scheme. 
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DELPHI cross-section limits 
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Figure 28: Expected cross-sections in pb at 208 GeV (dots) versus the chargino mass in a) 
for AM> 10 GeV/c^ and b) for AM~ 3 GeV/c^. The spread in the dots originates from 
the random scan over the parameters /i and M2. A heavy sneutrino (mj>> 1000 GeV/c^) 
has been assumed in a) and m,>>M-± in b). The minimal cross-sections below the mass 
limits are indicated by the horizontal lines. 
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DELPHI Xi mass limits 



^105 
o 

^02.5 

+ ^ 
^100 



97.5 - 



95 



92.5 




90 



87.5 - 



85 



82.5 - 



80 



Expected mass limit 
Obtained mass limit 



3 GeV/c^ 



10 



20 



30 



Kinematic limit 



Stable X? 

m : > 1000 GeV/c^ 



40 50 60 

Mxl-M%? (GeV/c') 



Figure 29: The chargino mass limit as function of the AM value under the assumption of 
a heavy sneutrino. The limit applies to the case of a stable Xv The straight horizontal 
line shows the kinematic limit. 
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Figure 30: The chargino mass limit as function of M2 for tan (3=1, under the assumption 
of a heavy sneutrino. The straight horizontal line shows the kinematic limit of the 
chargino production. The limit applies in the case of a stable Xv 
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Figure 31: Region excluded in the plane of the mass of the lightest neutralino versus that 
of the lightest chargino under the assumption of a heavy sneutrino, for tan (3 = 1.0, 1.5 
and 35. The thin lines show the kinematic limits in the production and the decay. The 
dotted line (partly hidden by the shading) shows the expected exclusion. The dashed 
region is not allowed in the MSSM. The limit applies in the case of a stable Xi- The mass 
limit on the lightest neutralino is indicated by the horizontal dashed line (for a discussion 
on the low mass region see section 6.3). 
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Figure 32: Efficiency for selecting a single almost stable chargino, as function of its mass, 
at the centre-of-mass energies of the years 1999 and 2000. 
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Figure 33: (a) Efficiency for selecting a single 75 GeV/c^ chargino in the search for 
displaced decay vertices (kinks) at the centre-of-mass energy of 196 GeV, as function of 
its decay radius, (b) Trigger efficiency for the selected charginos. (c) Trigger efficiency 
for all 75 GeV/c^ charginos, whether or not they were selected. 
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Figure 34: Selection times trigger efficiency in the search for nearly mass-degenerate 
charginos in the search with an ISR photon. The efficiency for higgsinos at a/s = 206 GeV 
at the different chargino masses and AM values fully simulated is given as an example. 
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Figure 35: Expected number of background events (top) and number of selected events in 
the real data (bottom) after the final selection, for the year 2000 data with the TPC fully 
operational, as function of the points in the plane (M^±, AM) considered in the chargino 
search with an ISR photon. 
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Figure 36: Regions in the plane (M-+,AM) excluded using: the search for high AM 
charginos; the search for soft particles accompanied by ISR; and the search for long-lived 
charginos. The two scenarios are (upper plot) the one in which the lightest chargino is 
a higgsino and (lower plot) the one in which the lightest chargino is a gaugino. For the 
second scenario, the limits are valid in the heavy sfermion approximation, while for the 
higgsino scenario it is sufficient that M? > M-+. 
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Figure 37: Contour plots of the upper limits obtained on the cross-sections for X1X2 
production at ^/s = 206 GeV. The data at all energies were used, assuming the 
cross-section energy dependence expected at a chosen high mo point in the higgsino region 
where the neutralino searches play an important role (mo=l TeV/c^, /i=-60 GeV/c^, 
M2=200 GeV/c^). In each plot, the different shadings correspond to regions where the 
cross-section limit in picobarns is below the indicated number. For figures a), b), c), X2 
decays into x? and a) e+e~, b) and c) qq, while in d) the branching ratios of the 

Z were assumed, including invisible states. The dotted lines indicate the kinematic limit 
and the defining relation M^o > M^o. 
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Figure 38: Upper limits on the cross-sections for X2Xi production with Xi'~^X2lQ (^=354) 
including data up to ^/s =208 GeV. The different shadings correspond to regions where 
the cross-section limit in picobarns is below the indicated number. The X2 "^^s assumed 
to decay 100% into XiQQ (Isft plot), and into Xi7 (right plot). The limits in the left 
plot are based on the acoplanar jets and multijets selections, while those in the right plot 
derive from the search for multijets with photons. 
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Figure 39: Excluded regions in the (/U,M2) plane for tan/5 = 1 for mo = 1000 GeV/c^. 
The shaded areas show regions excluded by searches for charginos and the hatched areas 
show regions excluded by searches for neutralinos. The thick dashed curve shows the 
isomass contour for M^o = 39.2 GeV/c^, the lower limit on the LSP mass obtained at 
tan (3=1. The chargino exclusion is close to the isomass contour for M-± at the kinematic 
limit. From chargino searches alone the lower limit on M^o is M^o = 38.2 GeV/c^. 
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Figure 40: The lower limit on the mass of the lightest neutralino, Xi, as a function of tan P 
assuming a stable Xv The solid curve shows the limit obtained for mo =1000 GeV/c^, 
the dashed curve shows the limit obtained allowing for any mo assuming that there is no 
mixing in the third family {Ar = /itan P, Af, = /itan /3, At = yu/tan P), and the dash-dotted 
curve shows the limit obtained for any ttiq allowing for mixing with AT-=Af,=At=0. The 
steep solid (dashed) line shows the effect of the searches for the Higgs boson for the 
maximal M-^o scenario (no mixing scenario), mo < 1000 GeV/c^ and Mt= 174.3 GeV/c^, 
which amounts to excluding the region of tan/3 < 2.36(9.7). 
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Figure 41: The masses as a function of tan/5 of the hghtest stau (sohd curve), stop 
(dotted curve) and sbottom (dashed curve), at the largest allowed \fj,\ for the smallest 
non-excluded M2 value. Mass splitting in the stau (sbottom, stop) sector in the form 
Ar — /itan/3 ( Af, — yUtan/3, At — /i/tan/3) was assumed, with Ar = {Af, = At =0). 
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Figure 42: The limit set on the LSP mass by chargino (dotted hne) and neutrahno 
searches (dashed hne) as a function of the sneutrino mass in the case when the hghtest 
stau is degenerate in mass with the hghtest neutrahno. The combined hmit is shown as 
sohd hne. These hmits are vahd for any model of stau mixing. In any particular model, 
in which mixing in the stau sector and in the sbottom and stop sector is related a more 
stringent limit on the LSP mass can be set. 
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Figure 43: The minimum sneutrino mass (thin dashed curve) and the e^ mass (thin 
dotted curve and hght shading) allowed by the slepton and neutralino searches, as a 
function of tan together with the limits on the chargino mass (thick solid curve and 
dark shading, and thick dash-dotted curve). The chargino mass limit indicated by the 
solid curve and the sneutrino and selectron mass limits were obtained assuming no mass 
splitting in the third sfermion family [A^. — /itan/5=0 in particular). The selectron mass 
limit is valid for Mgj, — M^o > 10 GeV/c^. The chargino mass limit indicated with the 
dash-dotted curve was obtained allowing for mass splitting in the third sfermion family, 
with Ar = Ab = A=0. 
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Figure 44: Regions excluded at in the (/i, M2) plane for tan/3 = 35 and different as- 
sumptions of mixing in the third family. Light shaded regions are excluded by searches 
at LEPl energies, darker shading bounded by the solid line marks regions which are 
excluded by neutralino searches, while intermediate shaded regions bounded by the thin 
dashed line are excluded by chargino searches. Mixing terms in the form — /ztan f3 
(Afc — yutan/?. At — fi/tanP) were considered. For plot a) Af, = At = A^- = was assumed. 
The no-mixing scenario was used in b) {A^ — fitan P=0, Aj, — fitan P=0, — /i/tan/3=0). 
Plot a) is for the mo values giving the lowest non-excluded LSP mass. In the vertically 
(horizontally) hatched areas the stau (the sbottom) is the LSP. Relevant isomass contours 
of the lightest neutralino are also shown (M^o=46 GeV/c^ in a) and M^o=50 GeV/c^ in 
b)). 



